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ABSTRACT 


)  A  new  geometric  adaptive  control  (GAC)  system  has  been  developed 
for  Improving  the  accuracy  and  stability  In  machining  cylindrical 
workpiece.  The  system  takes  Into  consideration  error  generation  process, 
machine  tool  dynamics,  and  metrology.  The  new  stochastic  model 
constructed  In  this  thesis  provides  not  only  Insight  Into  the  nature  of 
machining  processes,  but  also  a  means  for  the  systematic  design  of 
geometric  adaptive  controllers. 

The  adaptive  control  algorithms  for  the  GAC  system  are  derived  by 
using  the  theory  of  self-tuning  control  (STC).  The  algorithms  are  simple, 
robust,  and  suitable  for  microprocessor  Implementation.  Whltehouse's 
multiprobe  measurement  Is  modified  for  use  In  the  control  system  for  the 
first  time.  Interesting  properties  of  this  measurement  In  control  are 
analyzed  and  predicted. 

The  results  of  simulation  show  that  the  GAC  system  can  Improve 
both  the  accuracy  and  the  stability  considerably.  Through  the  theoretical 
analysis  It  was  possible  to  resolve  the  control  problems  Imposed  by 
multiprobe  measurement,  and  to  achieve  near  optimum  control 
performance.  A  comparison  with  another  GAC  system  which  uses  FCC 
technique  Is  made  to  show  that  the  proposed  GAC  system  does  have 
better  performance  than  other  existing  GAC  systems  In  precision 


machining. 


CHAPTER  1 


INTRODUCTION 


The  progress  of  machining  accuracy  has  been  accelerated  by  the 
requirement  of  high  precision  parts  in  space  exploration,  electronics 
miniaturization,  and  laser  applications  [70,93].  Not  only  do  these 
technologies  need  precision  hardware,  which  are  usually  produced  by 
machining,  to  accomplish  their  extraordinary  performance.  In  the 
traditional  industries  the  pursuit  for  better  accuracy  never  stops. 
Accurately  machined  parts  can  assure  assembly  efficiency  in  mass 
production  systems,  reduce  noise  and  vibration  for  improving  the 
quality  of  the  product,  and  simplify  product  design  [18]. 

It  can  be  seen  that  cylindrical  machined  parts  exist  in  most,  if 
not  all,  of  the  products.  This  is  particularly  true  for  those  having 
rotary  motions.  Dimensional  accuracy,  roundness,  and  finish  are  the 
three  most  important  aspects  of  the  accuracy  of  cylindrical  parts. 
Dimensional  accuracy  refers  to  the  bias  error  of  the  effective 
diameter  of  the  machined  part  from  the  specified  diameter  in  the  part 
drawing.  It  determines  the  clearance  or  tightness  in  the  assembly  of 
a  hole  and  a  shaft.  Dimensional  accuracy  is  usually  specified  in  the 
design  stage  as  tolerance  to  allow  for  the  machining  errors. 

It  has  been  shown  by  Moore  [70]  that  roundness  becomes 
proportionally  more  critical  as  tolerances  become  tighter. 


Roundness,  or  more  precisely  out-o f -roundness ,  refers  to  the 
deviation  of  the  cross-section  of  a  cylindrical  workpiece  from  an 
ideal  circle.  Depending  on  the  types  of  applications,  out-of- 
roundness  may  result  in  various  kinds  of  problems,  such  as  leakage, 
stress  concentration,  noise  and  vibration,  and  so  on  [34].  host  of 
the  time,  out-of-roundness  is  regarded  as  the  form  error  of  low 
frequency  undulations  on  the  part  surface. 

Unlike  dimensional  accuracy  and  roundness,  finish  is  a 
microscopic  description  of  the  part  surface.  The  machined  part 
usually  has  many  random  and  high  frequency  undulations  on  the  surface 
in  addition  to  out-of-roundness.  The  qualitative  interpretation  of 
these  small  undulations  through  human  sensors,  vision  and  touch,  is 
finish.  Its  effects  on  the  performance  of  the  product  are  also  of 
microscopic  nature,  such  as  friction,  wear,  erosion,  and 
ref lectability. 

According  to  the  above  descriptions,  it  is  understandable  that 
the  development  in  cylindrical  machining  has  always  been  the  most 
rewarding  in  the  history  of  machine  tools.  Lathes  and  drilling 
machines  are  used  for  machining  basic  cylindrical  surfaces.  Other 
machine  tools  such  as  boring,  reaming,  grinding,  and  honing  are 
devised  to  trim  excessive  dimensions,  to  true  round  shapes,  or  to 
polish  surfaces.  However,  the  higher  the  precision  the  process  can 
achieve,  the  higher  the  unpredictability  it  has,  and  the  higher  the 
skill  it  needs  [53].  Some  automatic  control  systems  for  geometrical 
accuracy  are  needed  to  overcome  this  barrier.  for  setting  up 
effective  control  systems,  it  would  be  helpful  to  identify  the 


sources  of  and  the  processes  generating  geometrical  inaccuracies  or 

errors  on  the  machined  part  in  advance. 

1.1  Sources  of  Geometrical  Errors 

The  sources  of  geometrical  errors  can  be  categorized  as  follows: 

i.  Structural  inaccuracies  and  deformation  -  All  members  of  the 
machine  tool  are  produced  by  using  various  production 
processes.  Consequently  they  are  subject  to  manufacturing 
errors.  And  because  of  the  finite  rigidity  of  the  members, 

structural  deformation  is  induced  whenever  there  are  forces  of 
cutting,  clamping,  gravity,  or  heat  expansion.  These  errors  are 
copied  to  the  machined  surfaces  through  tne  machining  processes. 

ii.  Vibrations  -  Vibration  is  the  most  annoying  problem  for  the 
machine  tool  builders  and  operators.  It  can  be  from  anywhere, 
unbalanced  rotating  parts,  loose  objects,  or  the  floor  of 
machine  base.  The  most  detrimental  vibration  is  from  the 
cutting  process  itself,  called  chatter.  Vibration  not  only 
reduces  the  accuracies  of  the  machined  part  but  also  shortens 
the  lives  of  the  tool  and  the  machine. 

iii.  Machining  processes  -  Besides  the  process  related  chatter,  the 
tool  may  leave  unwanted  marks  on  the  machined  surfaces  as  it 
moves  relative  to  the  workpiece.  Also  the  wear  of  the  tool 
affects  finish  and  dimensional  accuracy. 

iv.  Microstructural  changes  -  This  involves  the  rupture  or  the 
plastic  deformation  of  the  metallic  grains  of  workpiece  caused 


by  the  cutting  forces  and  the  effect  of  friction  on  the  tool- 
workpiece  interface. 

Roughly  speaking,  the  first  two  categories  contribute  to  dimensional 
error  and  out-of-roundness .  And  the  last  two  categories  have  major 
effects  on  surface  finish.  It  is  not  wise  to  consider  all  the  error 
sources  in  developing  an  automatic  control  system  for  geometrical 
accuracy.  Some  error  sources  are  of  minor  importance,  and  can  be 
neglected.  Some  of  them  may  never  be  good  for  automatic 
compensation.  Therefore,  it  is  a  good  practice  to  learn  about  the 
available  techniques  of  error  reduction.  Then  we  can  determine  the 
error  sources  which  should  be  considered  in  the  control  system. 

1.2  Methods  for  Error  Reduction 

According  to  Blaedel  [15],  the  methods  for  error  reduction  can 
be  divided  into: 

(a)  error  avoidance 
*eliminating  the  error  sources 
*altering  the  processes  of  error  generation 

(b)  error  compensation 
*precalibrated 

*active 

The  methods  of  error  avoidance  prevail  in  most  of  the  machine  tool 
builders.  They  typically  need  painstaking  trials  and  tests  in  order 
to  identify  the  true  error  sources  or  error  generation  processes. 
Then,  one  can  reduce  errors  by  redesigning,  modifying,  or  carefully 
manufacturing  the  critical  members  of  the  machine  tool.  For  cutting 
process  related  errors,  sometimes,  it  is  necessary  to  limit  cutting 


conditions  at  the  expense  of  productivity.  In  general,  error 
avoidance  applies  to  those  error  sources  which  are  steady  or  which 
can  be  eradicated  once  for  all. 

Error  compensation,  on  the  other  hand,  is  suitable  for  dynamic 
errors.  It  attempts  to  model  the  process  of  error  generation  from 
the  measured  variables  and  then  to  predict  the  control  inputs  for 
compensation.  If  the  measurement  is  performed  before  or  after  the 
cutting  operations  and  then  the  errors  are  compensated  for,  it  is 
called  precalibrated  error  compensation.  The  principle  of 
precalibrated  compensation  is  based  on  the  existence  of  systematic 
errors  for  modeling.  The  off-line  measurement  simplifies  the  system 
of  precalibrated  compensation.  However,  certain  nonrepeatable  or 
dynamic  errors,  chatter,  for  instance,  will  never  be  sensed  and 
corrected.  Production  time  is  also  wasted  in  measurement  for 
calibration. 

Active  error  compensation,  with  in-process  measurement, 
therefore,  emerges  as  the  most  promising  method  for  reducing 
unpredictable  errors  [47,48,77,78].  One  type  of  in-process 
measurement  is  to  measure  the  error  motions  of  certain  critical 
members  of  the  machine  by  attaching  external  masters  to  those 
members.  The  other  type  is  to  measure  the  geometrical  profile  of  the 
machined  surface  directly.  Although  sometimes  the  latter  type  of 
active  compensation  has  accessibility  problems  in  setting  up  the 
instruments,  its  performance  is  more  effective  than  that  of  the 
former  type  because  of  the  direct  assessment  of  geometrical  errors. 
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The  modeling  of  Che  process  in  active  compensation ,  as  shown  in 
Fig.  1.1,  is  quite  complex.  The  unpredictable  properties  of  the 
system  can  be  resolved  by  devising  an  adaptive  model  for  the 
compensator  or  controller  according  to  in-process  measurement , 
compensation  input,  and  appropriate  control  criteria.  This  has  been 
termed  "geometrical  adaptive  control"  (GAC)  [77,48,47]. 

1.3  GAC  System  and  Related  Problems 

Basically,  a  GAC  system  should  consist  of  in-process 
measurement,  on-line  modeling  or  system  identification,  and  real  time 
control.  It  is  best  for  the  applications  in  fine  turning,  honing, 
and  cylindrical  grinding  processes.  Like  its  companions  in 
machining,  adaptive  control  for  optimization  (ACO)  and  adaptive 
control  for  constraint  (ACC)  [99] ,  the  development  of  GAC  has  been 
relatively  slow.  Until  now,  only  a  handful  of  recent  research  papers 
can  be  seen  [50,84,55,68,69,87,102].  Most  of  these  research  work 
still  fail  to  fully  resolve  the  following  problems: 

A.  Problems  in  measurement  -  Apart  from  the  problems  of  the 
accessibility  and  reliability  of  sensors  in  harsh  production 
environments,  a  major  concern  in  GAC  is  the  separation 
problem.  That  is,  with  the  workpiece  attached  to  the  rotating 
spindle,  any  in-process  measurement  of  the  surface  geometry  will 
be  contaminated  by  the  spindle  error  motions,  which  are 
difficult  to  be  separated.  The  external  masters  used  in  [84,55] 
can  measure  and  compensate  for  the  spindle  error  motion  only. 
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Figure  1.1  Geometric  Adaptive  Control  System 


A.  Problems  in  system  '  modeling  -  A  thorough  understanding  of 
machining  processes  has  been  elusive  for  the  inconsistencies  of 
experimental  findings  caused  by  different  cutting  conditions, 
tool  geometries,  and  work  materials.  Traditional  models  are 
deficient  in  handling  the  nonlinear,  nonstationary,  and 
stochastic  nature  of  machining  processes.  On  the  other  hand, 
the  stochastic  time  series  model  used  by  Wu  and  his  coworkers 
(84,55,68,69]  bears  little  physical  meaning. 

C.  Problems  in  adaptive  control  -  In  the  absence  of  appropriate 
models  for  machining  processes ,  most  of  the  adaptive  control 
techniques  of  GAC  are  quite  crude  and  inadequate  in  assuring  the 
stability  of  the  GAC  system.  A  systematic  approach  for  the 
design  of  adaptive  controller  is  still  lacking. 

1.4  The  Objective  of  the  Research 

The  objective  of  this  thesis  is  to  propose  a  new  GAC  system 
which  can  improve  the  precision  of  today's  cylindrical  machining 
processes  by  at  least  an  order  of  magnitude.  To  realize  this 
objective  the  above  problems  should  be  tackled  more  rigorously. 

In  Chapter  2,  we  Investigate  in  detail  the  nature  of  the 
problems  in  measurement  and  modeling.  The  basic  principles  and 
characteristics  of  metrology  are  studied  to  help  develop  a  suitable 
technique  of  in-process  measurement.  For  the  purpose  of  modeling,  we 
review  machine  tool  dynamics  and  try  to  find  the  most  significant 
nature  of  machining  processes  for  developing  a  generic  but  physically 
meaningful  model  of  the  GAC  system. 


In  Chapter  3,  the  general  theory  of  self-tuning  control  is 
introduced.  It'  is  applicable  to  our  GAC  system  for  its  ability  to 
control  time  varying  and  stochastic  systems.  Also,  it  satisfies  the 
requirement  for  a  systematic  design  of  adaptive  controllers.  In 
general,  self-tuning  control  consists  of  two  parts,  identification 
and  control.  There  are  many  algorithms  for  parameter  identification 
and  controller  design.  To  develop  an  effective  self-tuning 
controller,  some  practical  1  sues  have  to  be  considered  cautiously  in 
the  construction  of  algorithms. 

Mathematical  models  of  the  GAC  system  and  the  corresponding 
adaptive  control  algorithms  in  discrete  forms  for  digital  computation 
are  presented  in  Chapter  4.  The  derivations  are  based  on  the 
investigations  in  Chapters  2  and  3.  The  stability  and  convergence 
properties  of  the  control  algorithms  are  analyzed.  In  addition, 
special  properties  of  multiprobe  measurement,  a  method  of  in-process 
measurement,  in  the  GAC  system  are  studied. 

To  evaluate  the  performance  of  the  developed  algorithms  in 
Chapter  4,  simulation  is  carried  out  and  the  results  are  presented 
and  discussed  in  Chapter  5.  The  determination  of  some  important 
factors  associated  with  the  practical  issues  in  self-tuning  control 
can  be  observed  through  the  simulation  results.  The  geometrical 
accuracies  and  the  increase  in  stability  are  the  measures  of  the 
performance  of  the  GAC  system.  The  simulation  also  accomplishes  the 
interesting  studies  of  multiprobe  measurement. 


In  Chapter  6  comparisons  with  Forecasting  Compensatory  Control 
(FCC),  which  has  been  used  for  roundness  control  [84,55],  are  given 
to  show  that  the  GAC  system  we  developed  here  does  have  significant 
improvement  in  performance,  accuracy  and  stability  over  other 
existing  techniques.  We  make  the  comparison  in  theoretical  analysis, 
digital  simulation,  and  analog  test. 

Chapter  7  concludes  the  work  we  have  done  and  suggests  future 
work  which  will  support  the  success  of  this  new  development. 
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CHAPTER  2 

METROLOGY  AND  DYNAMICS 

The  development  of  GAC  systems  involves  the  consideration  of 
metrology,  machine  tool  dynamics,  and  control  theory.  In  this 
chapter  the  discussion  will  be  focussed  on  metrology  and  dynamics. 
The  control  theory  will  be  presented  in  next  chapter. 

2.1  Metrology 

"Metrology  is  the  science  of  measurement,  which  determines  a 
dimension  through  the  use  of  some  type  of  calibrated  device  that 
permits  the  magnitude  of  the  dimension  to  be  determined,  directly  and 
indirectly,  in  scalar  units"  (31].  The  evolution  of  metrology  in 
manufacturing,  since  the  early  nineteenth  century,  has  gone  through 
three  stages:  passive,  active,  and  dynamic  [80].  In  the  first  stage, 
passive  metrology  only  served  to  reduce  assembly  work  and  to 
guarantee  interchangeability  by  checking  the  dimensional  parameters 
after  manufacturing.  In  the  next  stage,  active  metrology  was 
introduced  in  order  to  achieve  quality  and  interchangeability  more 
economically.  The  measurement  is  brought  nearer  to  the  manufacturing 
process  so  that  inaccuracies  can  be  actively  corrected  by  a  certain 
precalibrated  compensation  technique  as  mentioned  in  Chapter  1.  In 
the  third  stage,  starting  from  1960's,  dynamic  metrology  emerged  as 
an  essential  part  of  the  integrated  manufacturing  system  for  a 


continuous  and  in-process  measurement.  The  measurement  can  be  used 
not  only  to  correct  the  inaccuracies  but  also  to  optimize  the 
production  system  by  the  feedback  or  feedforward  control  action. 
When  correction  of  the  geometrical  inaccuracies  is  the  major 
function,  the  integrated  system  is  the  GAC  system  as  defined  in 
Chapter  1;  otherwise,  it  can  be  an  ACC  or  ACO  system. 

In  parallel  with  the  above  evolution,  the  instrument  for 
metrology  has  also  developed  from  James  Watt's  micrometer  to  today's 
highly  sophisticated  and  precise  laser  interferometers  [44]. 
Meanwhile,  the  advent  of  electronics  and  digital  computers  in 
metrology  have  made  the  measurement  more  consistent  and  efficient. 
In  addition,  metrology  also  extends  from  simple  distance  measurement 
to  surface  profile  assessment.  More  metrological  parameters  are 
needed  to  quantify  the  functions  or  performances  of  a  manufactured 
part.  Standards,  such  as  ANSI,  BS,  and  ISO,  have  various  documents 
on  metrology  [42].  They  are  used  to  justify  the  measurement  of  the 
metrological  parameters. 

Metrology  in  manufacturing,  in  general,  can  be  categorized  into 
three  types:  workpiece  metrology,  machine  tool  metrology,  and  in- 
process  metrology.  Although,  for  the  GAC  system,  we  are  interested 
in  in-process  metrology,  measuring  workpiece  during  machining 
operation,  the  basic  metrology  principles  and  instrumentation  of  the 
other  two  types  are  presented  here  to  help  in  identifying  and  solving 
problems  in  in-process  metrology  or  measurement.  Since  we  are 
proposing  a  GAC  system  for  cylindrical  machining  processes,  most  of 
the  discussion  will  be  confined  to  cylindrical  machined  parts  and  the 
associated  machine  tools. 


2.1.1  Workpiece  Metrology 


Workpiece  metrology  is  the  easiest  to  conduct,  as  the  workpiece 


can  be  placed  in  an  ideal  metrology  room  and  measured  carefully.  For 


a  cylindrical  workpiece,  we  are  most  interested  in  the  measurement  of 


roundness.  The  methods  of  measuring  roundness  can  be  divided  into 


production  type  and  absolute  type  [70].  The  methods  of  production 


type  are: 


(i)  diametral 


(ii)  circumferential  confining  gage 


(iii)  rotating  on  centers 


(iv)  V-block 


(v)  three-point  probe 


These  methods  are  relatively  inexpensive,  easy  to  operate,  and  rugged 


in  use.  The  disadvantage,  however,  is  that  the  relative  measurement 


of  roundness  can  only  be  used  for  passing  or  rejecting  the  part  in 


production  line.  The  method  of  rotating  on  centers  is  even 


susceptible  to  errors  from  both  the  measuring  setup  and  the  part 


itself. 


The  major  criticism  of  the  diametral  method  and  the  V-block 


method  is  their  inability  to  detect,  for  example,  the  iso-diametrical 


out-of-roundness  as  shown  in  Fig.  2.1.  From  the  configurations  in 


Fig.  2.2,  the  out-of-roundness  is  determined  by  the  magnification 


factor  [ 90] . 
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Figure  2.1  Isodiametrical  Shapes  [1] 


Figure  2.2  Configurations  of  V-block.  for  Measuring 
Ouc-of -Roundness 


where  a  is  the  angle  of  V-block  and  n  is  the  number  of  lobes 

distributed  on  the  part  surface  uniformly.  F^  equal  to  zero  means 

that  the  out-of -roundness  will  not  be  detected  for  an  n-lobed  part. 

It  can  be  shown,  from  (2.1),  that  the  diametral  gage,  equivalent  to  a 

180°  V-block,  is  unable  to  detect  the  out-of-roundness  of  the  part 

with  odd-numbered  lobes.  And  60  V-block  is  not  good  for  5,  7-lobed 

parts;  90°  V-block  not  for  7,  9-lobed  parts  [70].  More  specifically, 

the  V-block  of  Fig.  2.2(a)  can  not  detect  the  out-of-roundness  of  the 
* 

n  -lobed  part ,  where 

n*  -  2*k/8  ±  1  (2.2) 

is 

and  k  is  any  positive  integer  making  n  integer.  The  same  property 
can  al6o  be  found  in  the  multiprobe  measurement  of  Chapter  4. 

The  configuration  of  V-block  method  is  also  analogous  to  that  of 
centerless  grinding.  As  shown  in  Fig.  2.3,  the  edges  of  V-block 

correspond  to  the  grinding  wheel  and  the  control  wheel,  and  the 

indicator  above  the  V-block  corresponds  to  the  workplate  of  the 
grinding  machine.  Therefore,  (2.1)  or  (2.2)  explains  why  the  ground 
workpiece  usually  has  odd-numbered  lobes,  and  the  out-of-roundness 
diminishes  as  the  workplate  raises  to  a  higher  position. 

The  absolute  method  means  that  the  roundness  is  assessed 
directly  from  the  surface  profile  measured  in  a  roundness  measuring 
machine  which  has  an  extrinsic  datum  [34],  The  measuring  machine  can 
be  either  with  overhead  spindle  or  with  rotating  table  as  shown  in 

Fig.  2.4.  The  overhead  spindle  type,  such  as  Moore  Universal 


Measuring  Machine,  has  a  stylus  tracing  around  the  cylindrical  part 
on  a  fixed  flat  table.  It  is  good  for  heavy  and  short  parts.  The 
rotating  table  type,  such  as  Talyrond,  on  the  other  hand,  is  good  for 
light  parts.  The  stylus  assembly  positioned  on  a  fixed  stand  can  be 
adjusted  freely.  Specifications  on  stylus  tip  radius,  stylus  static 
force,  cycles  per  revolution,  and  other  operational  information  can 
be  found  in  British  Standard  [17]  and  ANSI  Standard  [1]. 

The  movement  of  the  stylus  is  amplified  and  plotted  on  a 
circular  chart.  The  two  most  popular  methods  for  evaluating  the 
out-of-roundness  are:  minimum  zone  center  (MZC)  and  least  squares 
center  (LSC)  [17].  The  MZC  method,  as  shown  in  Fig.  2.5,  needs 
several  graphical  trials  done  manually  to  determine  the  two 
concentric  circles  having  the  minimum  radial  difference.  This  radial 
difference  is  defined  as  out-of-roundness .  Recently,  Murthy  and 
Abdin  [72]  devised  the  simplex  search  for  finding  the  MZC  in  a  more 
systematic  way. 

The  LSC  method,  as  shown  in  Fig.  2.6,  computes  and  finds  the 
center  and  the  radius  of  the  least  squares  circle  by 


2 


2 


(2.3) 


The  out-of-roundness  is  the  summation  of  the  largest  outward 


Lust  squares  circle 


1 


deviation  of  the  polar  profile  from  the  LS  circle  and  the  largest 
inward  deviation.  Since  the  LS  circle  can  be  determined  uniquely, 
this  method  is  more  suitable  for  computer  controlled  measuring 
machines . 

Special  care  should  be  taken  in  assessing  out-of-roundness  when 
the  cylindrical  part  is  miscentered  with  respect  to  the  axi6  of 
rotation  of  the  measuring  machine.  Fig.  2.7  shows  that  a  perfect 
round  part  will  have  a  limacon-shaped  polar  graph  provided  that 
eccentricity  exists.  The  reason  for  this  kind  of  distortion  is  the 
suppressed  radius  used  in  plotting  the  polar  graph  [110].  The 
suppressed  radius  is  the  preset  distance  between  the  stylus  tip  and 
the  axis  of  rotation.  It  should  be  as  close  to  the  mean  diameter  of 
the  part  as  possible  so  that  the  undulation  of  the  part  surface  can 
be  distinguished  by  a  high  amplification  factor. 

Therefore,  the  above  descriptions  on  assessing  out-of-roundness 
are  just  approximately  true  if  the  eccentricity  ratio  is  less  than 
15%  [95,1].  Whitehouse  [110],  Chetwynd  [22],  and  McCool  [63] 
proposed  a  more  accurate  approach  for  numerical  assessment.  This 
approach  uses  the  following  regression  equation: 

r  ■  R  +  a  cos  6^  +  b  sin  6^  +  e^ 

i  -  1,2 . N  (2.4) 

By  tf i  least  squares  method,  the  estimated  values  for  R,  a,  and  b  are 
the  same  as  those  in  (2.3),  and 
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Figure  2.7  Umacon-shaped  Polar  Profile  of 
Mis-centered  Parc  [11 
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p(0^)  •  R  +  a  cos  9^  +  b  sin  0^ 

represents  the  limacon  graph  of  the  LS  circle  of  a  mis-centered  part. 
The  residual  error  e^  is  the  quantity  used  for  assessing  roundness. 
As  defined  in  the  British  Standard,  the  LSC  out-of -roundness  is  the 
maximum  peak-to-peak  error  of  e^.  Besides  the  improved  accuracy,  the 
limit  on  maximum  allowable  eccentricity  is  less  restrictive  for  this 
technique.  It  simplifies  the  centering  operation  during  measurement. 

Also,  (2. A)  shows  that  the  limacon  consists  of  the  constant  and 
the  first  harmonic  of  the  Fourier  series,  and  e^is  the  sum  of  all 
higher  harmonics.  Hence,  the  removal  of  limacon  from  the  polar  graph 
is  equivalent  to  using  a  high-pass  filter  in  signal  processing.  This 
technique  is  frequently  used  in  machine  tool  metrology. 

2.1.2  Machine  Tool  Metrology 

Basically  machine  tool  metrology  is  the  same  as  workpiece 
metrology  with  the  absolute  method,  except  that  the  part  to  be  sensed 
is  a  perfect  round  master  attached  to  the  spindle  of  the  machine 
tool.  Any  measured  movement  can  be  attributed  to  the  spindle  error 
motion,  which  may  generate  uneven  surface  profile  on  the  machined 
part  [20,86].  Since  the  position  of  the  sensor  is  where  the  tool 
resides  normally,  the  spindle  error  motion  in  the  radial  direction 
relative  to  the  sensor  is  particularly  important  in  producing  the 
geometrical  inaccuracies  of  dimensional  error,  out-of-roundness ,  and 


roughness . 


Unlike  the  absolute  method  of  workpiece  metrology,  the  setup  for 
machine  tool  metrology  is  determined  by  the  machining  processes.  For 
turning  and  most  grinding  operations,  the  setup  for  measurement 
involves  the  replacement  of  the  stationary  cutting  tool  by  a 
noncontact  sensor  and  the  rotating  workpiece  by  a  spherical  master, 
as  shown  in  Fig.  2.8(a).  For  the  boring  process  with  rotating  tool, 
Fig.  2.8(b)  shows  the  setup  with  two  perpendicular  probes  for 
measuring  the  error  motion  at  the  angular  position  of  the  tool.  The 
setup  in  Fig.  2.9,  an  alternative  for  the  setup  in  Fig.  2.8(b),  is 
more  convenient  for  using  single  probe  149].  Other  methods  related 
to  the  above  two  cases  can  be  seen  in  the  survey  by  Bryan  and 
Vanherck  [20]  and  Murthy  et  al.  [73]. 

Because  of  the  nonrepetitive  nature  of  the  spindle  error  motion, 
the  measurement  is  normally  made  of  readings  for  several  revolutions 
as  shown  in  Fig.  2.10;  it  is  defined  as  the  total  error  motion  polar 
plot.  The  mean  contour  of  the  total  error  motion  polar  plot  is 
termed  as  the  average  error  motion  polar  plot,  and  the  difference 
between  the  total  error  motion  polar  plot  and  the  average  error 
motion  polar  plot  is  called  the  random  error  motion  [101].  Using  the 
LSC  method,  the  eccentricity  can  be  located  in  the  average  error 
motion  polar  plot. 

Since  the  spindle  runs  normally  at  the  operating  speed  during 
measurement,  the  induced  once-per-revolution  unbalanced  motion  may 
mix  with  the  misalignment  of  the  master.  One  of  the  approaches  for 
removing  the  misalignment  effect  from  the  measurement  is  to  take  a 
second  measurement  at  a  fairly  low  spindle  speed  to  get  a  new  average 
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8  Techniques  of  Machine  Tool  Metrology 

(a)  Stationary  Tool  Type 

(b)  Rotary  Tool  Type 


Another  Measurement  Technique  for  Process  of 
Rotating  Tool  Type  [49] 


Al  total  error  MOTION 


Figure  2.10  Typical  Measurement  of  Spindle  Radial 
Error  Motion  [101] 


motion,  then  to  find  the  eccentricity  purely  caused  by  misalignment. 
The  digital  technique  we  mentioned  is  good  for  this  type  of  signal 
processing. 

2.1.3  In-process  Metrology 

For  very  precise  measurement,  we  can  no  longer  assume  perfect 
spindle  in  workpiece  metrology  or  perfect  master  in  machine  tool 
metrology.  This  situation  corresponds  to  in-process  metrology  which 
involves  the  measurements  of  both  the  wortcpiece  geometry  and  the 
spindle  error  motion.  The  reversal  method  [32,1],  as  shown  in  Fig. 
2.11,  separates  the  spindle  error  motion  e(0)  and  the  surface  profile 
s(0 )  by 


e(0)  -  CvL(0)  -  v2(0))/2 
s(8)  =  (v^©)  +  v2(0))/2 

where  v^(0)  is  the  measurement  taken  at  position  1  and  v2(0)  at 
position  2.  Whitehouse  [111]  devised  the  multi-step  method  as  shown 
in  Fig.  2.12.  In  step  i,  the  measurement  is 

vi(0)  «  s(0  +  a±)  +  e( 0 ) ,  i  -  0,1,... ,n. 

Through  appropriate  linear  combination 


c(0)  *  Ea.v, (0)  “  Ea.s(0  +  a.) 


Figure  2.11  Reversal  Method  for  In-Process  Measurement 


Figure  2.12  Multi-Step  Method  for  In-Process  Measurement 


s(0)  can  be  evaluated  from  the  Fourier  expansion  of  c(6).  Although 
the  reversal  method  has  more  mechanical  problems,  the  multi-step 
method  suffers  harmonic  distortion.  Both  of  these  methods  assume 
that  the  spindle  error  motion  is  repetitive,  which  is  not  generally 
true  in  practice.  Thus,  Whitehouse  [111]  proposed  the  multiprobe 
method  for  non-repetitive  spindle  errors.  This  is  particularly  good 
for  our  GAC  system.  It  will  be  discussed  in  Chapter  4. 

2.2  Machine  Tool  Dynamics 

As  discussed  in  Chapter  1,  there  are  many  error  sources  which 
contribute  to  the  inaccuracies  of  a  machined  part.  Machine  tool 
dynamics  plays  a  very  important  role  in  transforming  the  error 
sources  into  the  geometrical  inaccuracies.  Also,  in  developing  a  GAC 
system  for  reducing  the  inaccuracies,  the  dynamic  stability  of  the 
machining  system  should  not  be  sacrificed,  however.  Thus,  it  is 
imperative  to  understand  machine  tool  dynamics  before  the  development 
of  the  GAC  system. 

In  general,  machine  tool  dynamics  deals  with  the  interaction 
between  the  cutting  process  and  the  machine  tool  structure  as  shown 
in  Fig.  2.13.  The  interest  in  the  investigation  of  machine  tool 
dynamics  can  be  attributed  to  the  most  obscure  and  delicate  problem 
in  metal  cutting  -  chatter  [94].  Chatter,  usually,  is  an  unwanted 
vibration  phenomenon  during  the  machining  operation.  It  spoils  the 
finish  and  the  dimensional  accuracy  of  the  workpiece,  shortens  the 
lives  of  the  tool  and  the  machine,  and  generates  uncomfortable  noise 


and  vibration. 


CUTTING 

PROCESS 
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There  are  two  basic  types  of  chatter:  forced  chatter  and  self- 
excited  chatter.  Forced  chatter  typically  results  from  periodic 
disturbances,  such  as  unbalanced  rotating  members,  intermittent 
cutting,  gear  backlash,  and  motions  transmitted  through  the  floor 
[16],  This  type  of  chatter,  though  important  in  practice,  relates  to 
fairly  simple  structural  dynamics  and  offers  no  insight  into  the 
dynamics  of  the  cutting  process.  Self-excited  chatter,  on  the  other 
hand,  occurs  only  when  there  are  any  machining  operations.  The 
interaction  of  cutting  dynamics  and  structural  dynamics,  under 
certain  conditions,  may  result  in  dynamic  instability  of  the 
machining  process.  Normally,  the  resulting  vibration  is  much  more 
violent  than  that  resulting  from  the  former  type.  Thus,  without  any 
confusion,  the  term  "chatter"  specifically  refers  to  the  self-excited 
type  in  the  following  text. 

2.2.1  Chatter  and  Cutting  Dynamics 

The  theory  of  chatter  in  metal  cutting  was  first  explored  by 
Arnold  [4],  He  proposed  the  falling  characteristic  of  force-speed 
relationship  for  explaining  the  chatter  phenomenon.  Although  this 
concept  is  still  quite  controversial  [40,98],  his  work  shed  some 
light  on  the  research  in  cutting  dynamics.  Basically,  the  cutting 
process  can  be  visualized  from  the  three  boundaries  in  Fig.  2.13(b) 
[30],  On  the  boundary  of  the  shear  plane,  the  length  and  the  angle 
of  the  plane  determine  major  cutting  forces  exerted  on  the  tool  faces 
[64],  On  the  boundary  of  the  tool-chip  interface,  the  tool  tip,  tool 
holder,  and  workpiece  may  deflect  owing  to  the  cutting  forces  on  the 
tool  rake  face.  The  geometry  and  the  tip  position  of  the  tool  are 


thus  varied.  Besides,  the  sliding  motion  of  the  chip  over  the  tool 
face  is  sometimes  complicated  by  the  built-up-edge  and  the  wear  on 
the  rake  face,  which  change  the  tool  geometry  further. 

Similar  observations  occur  on  the  boundary  of  the  tool-workpiece 
interface.  The  flank  wear  resulting  from  the  sliding  motion  changes 
the  effective  clearance  angle  of  the  tool  and  affects  the  reacting 
forces  on  the  boundary.  The  tool  and  the  workpiece,  again,  are 
deflected  by  the  reacting  forces.  Conceivably,  the  whole  machining 
process  is  quite  nonlinear  in  nature;  the  associated  factors  are 
Interrelated.  These  explain  the  frequent  inconsistencies  in  the 
experimental  findings  among  the  researchers. 

The  mechanisms  of  chatter,  generally,  can  be  divided  into  the 
velocity-dependent  type,  regenerative  type,  and  mode-coupling  type 
[104],  Although  they  act  simultaneously  in  practice,  regenerative 
chatter  is  recognized  to  dominate  most  of  the  unstable  vibrations  in 
metal  cutting  [43,52,97].  Furthermore,  Tobias  [98]  elaborated  the 
theory  of  regenerative  chatter  and  derived  a  stability  chart  which 
matches  very  well  with  the  test  results  as  shown  in  Fig.  2.14.  The 
lobing  characteristic  of  the  stability  chart,  which  is  common  to  most 
of  the  machining  processes,  can  hardly  be  explained  by  other  chatter 
theories.  These  facts  justify  the  regenerative  process  or  effect  as 
being  the  key  phenomenon  of  the  machining  process. 

Regenerative  effect  is  the  chip-thickness  variation  effect 
caused  by  the  surface  undulation,  y(t-t),  which  is  formed  in  a 
previous  cut  [98].  Hence,  as  seen  in  Fig.  2.13(b),  the  uncut  chip 


Vv 


thickness  for  a  cylindrical  machining  process  is  given  by 

6(t)  “  y(t)  -  py(t-x) 

where  x  is  the  time  required  for  each  revolution  of  the  spindle  and 
P,  overlap  factor,  is  introduced  to  account  for  the  degree  of  overlap 
in  successive  cuts,  p  is  equal  to  one  for  plunge  cutting,  but  zero 
for  thread  cutting.  In  general, 

0  <  p  <  i 

The  analytical  prediction  of  p  is  very  difficult;  nevertheless,  some 
attempts  have  been  made  by  Srinivasan  and  Nachtigal  [91]. 

Machine  tool  dynamics  can  be  conveniently  represented  by  the 
block  diagram  of  Fig.  2.15,  where  regenerative  effect  is  separated 
from  cutting  dynamics  for  its  particular  significance.  We 
deliberately  generalize  the  complexity  of  cutting  dynamics  by 
expressing  it  as  k^G^' .  k£  represents  the  static  cutting  stiffness, 
which  can  be  measured  in  steady  state  cutting.  The  dynamic  cutting 
stiffness  G  '  can  be  as  complex  as  the  theoretical  derivations  of 
Albrecht  [2]  and  Wu  [114] .  However,  linear  dynamics  as  expressed  by 

G  '(s)  -  1  +  Ks 
c 

is  adequate  to  account  for  the  physics  of  small  vibration  and  to 
accommodate  the  important  effect  of  penetration  rate  [98,89]. 
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As  for  structural  dynamics,  k  is  the  static  structural 

m 

stiffness  of  the  tool-workpiece-machine  system.  The  dynamic 

structural  stiffness,  G  ' ,  basically  has  the  transfer  function 

m 


G'(s) 
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which  is  a  multi-degree  of  freedom  system  with  damping  ratios  c^'s 

and  natural  frequencies  u)  ,'s.  If  there  is  any  mode-coupling  effect, 

ni 

G  '  can  be  expanded  into  a  transfer  matrix  to  correlate  the  vectors 
m 

of  cutting  forces  and  structural  deflections. 


A  further  simplification  of  the  representation  can  be  made  by 
letting 


G  -  G  'G  ' 
m  cm 

The  stability  borderline  in  the  stability  chart  can  be  constructed  by 
solving  the  characteristic  equation  of  the  system  in  Fig.  2.15 


1  +  _£  (1  -  Me  SX)Gm(s)  -  0,  s  -  ju>,  j  -  4-1 

m 

Thus 


k. 

,  .  m  1-p  cos  wt 

gR(ui)  “  “k - 2 - 

c  1  +  m  -  2ycos  <dt 


.  .  m  u  sin  u)i 

gl(u°  "  ~ - 2 - 

c  1  +  y  -  2p  cos  arc 


(2.5) 


(2.6) 
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where 


G  O)  -  g0(<*>)  +  jgT(<*>) 

OKI 

If  G  (jui)  and  u  are  known  in  advance,  the  chatter  frequency  w  can  be 
m  c 

found,  for  each  given  ,  by  solving 


p  sin  u  t 
_ c 

1  -  P  COS  U)  T 

c 

Then,  k  /k  can  be  obtained  from  either  (2.5)  or  (2.6). 
cm 

Fig.  2.16  is  the  stability  chart  for  the  special  case  when 


can  be  used 


4 


Figure  2.16  A  Theoretical  Stability  Chart  of  Fig.  2.15 

U*1 ,  u>  -  chatter  frequency,  w  -  natural 
frequency,  ft-  spindle  angular  speed 
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as  a  measure  of  the  stability  of  machining.  For  the  above  simple 

system,  with  typical  damping  ratio  of  C  *  0.05,  the  maximum  degree  of 

stability  is  k  /k  ■  0.105. 
c  m 

2.2.2  Methods  for  Chatter  Reduction 

There  are  many  ways  of  increasing  the  stability  of  machining  or 
chatter  control  [16,98].  Increasing  the  damping  in  the  vibration 
direction  is  the  most  obvious  one.  Increasing  the  structural 
stiffness  by  preloading  the  spindle  can  also  enhance  the  stability. 
In  the  shop  floor,  changing  the  cutting  conditions  is  is  frequently 
used  for  avoiding  chatter.  Small  depth  of  cut  or  width  of  cut,  which 
reduces  cutting  stiffness,  is  beneficial  in  stabilizing  the  machining 
process.  It  can  be  seen  from  Figs.  2.14  and  2.16  that  the  dynamic 
stability,  sometime,  can  be  achieved  by  using  high  cutting  speed 
which  exceeds  a  certain  critical  value.  Week  et  al.  [103]  even  used 
the  lobing  characteristic  of  the  stability  chart  to  avoid  chatter 
vibration  by  shifting  the  spindle  speed  to  the  stable  region. 

Most  of  the  above  methods  usually  restrict  the  productivity  or 
need  redesign  of  the  machine.  To  overcome  these  deficiencies,  active 
chatter  control  was  introduced  in  the  late  1960's.  The  system  of 
active  chatter  control  normally  has  a  feedback  controller  which  is 
designed  according  to  the  machine  tool  dynamics  as  shown  in  Fig. 
2.15.  The  feedback  system  of  Comstock  et  al.  [29]  measures  the 
relative  displacement  between  the  tool  and  the  spindle,  and  controls 
the  infeed  motion  of  the  tool  during  plunge  cutting.  Nachtigal  and 
his  coworkers  [74,75,60]  use  force  feedback  to  control  the  same 
machining  process. 


Later,  adaptive  active  chatter  control,  using  both  force  and 
displacement  feedbacks,  was  proposed  to  handle  the  problems  of  the 
variation  of  machine  tool  dynamics  during  operation  [76,67]. 
Basically,  these  methods  use  a  deterministic  model  of  machine  tool 
dynamics  for  stabilizing  the  machining  operation.  They  may  also 
attenuate  disturbances  through  the  increase  of  the  controller  gain 
[66].  However,  for  geometric  adaptive  control,  where  accuracy  is 
more  important  than  stability,  the  gain  increase  may  not  be  adequate 
for  the  unknown  nature  of  random  disturbances.  Therefore,  it  is 
essential  to  develop  a  stochastic  model,  based  on  the  similar 
reasoning  of  machine  tool  dynamics ,  to  control  the  random 
disturbances  more  effectively.  Also,  the  control  algorithms  must  be 
more  adaptive,  versatile,  and  easy  to  implement  for  automated 


machining  systems 


CHAPTER  3 


SELF-TUNING  CONTROL  THEORY 

Although  we  are  able  to  recognize  and  model  the  fundamental 
properties  of  machine  tool  dynamics  as  shown  in  Chapter  2,  the 
characteristic  parameters  of  the  system  can  never  be  assessed  exactly 
through  off-line  testings  and  analyses.  The  value  of  the  parameters 
may  change  with  the  type  of  working  material,  part  geometry,  tool 
geometry,  and  cutting  conditions.  Also,  the  changes  may  occur  during 
the  operation  because  of  cutting  force,  wear,  heat,  and  other 
disturbances.  Thus,  in  GAC,  we  need  an  adaptive  control  technique 
which  can  assess  the  system  characteristics  on  line  and  maximize  the 
performance  of  the  machining  process. 

As  shown  in  Fig.  1.1,  the  adaptive  controller  basically  consists 
of  two  important  functions: 

*on-line  identification  of  characteristic  parameters 
*real  time  control  for  optimizing  the  performance  index. 

An  ideal  or  optimal  adaptive  controller,  also  called  dual  controller, 
is  very  difficult  to  achieve  and  implement  due  to  the  nonlinear 
nature  of  the  interaction  between  the  identification  function  and  the 
control  function  [35,9]  The  interaction  is  nonlinear  in  the  sense 
that  the  controller  has  to  consider  the  future  uncertainties  of  the 
parameters  in  giving  the  control  inputs.  And  the  uncertainties  or 


probability  distributions  of  the  parameters,  called  hyperstate,  are 
functions  of  the  control  inputs  [112,5]. 

3.1  Self-tuning  Control  (STC) 

Self-tuning  control  is  a  suboptimal  adaptive  control  technique 
which  uses  the  certainty  equivalence  principle.  The  optimal 
controller  is  designed  by  considering  the  estimated  parameter  values 
as  the  true  parameter  values,  that  is,  by  ignoring  any  uncertainties 
of  the  estimated  parameters.  This  simplification  makes  the  adaptive 
control  system  linear  and  reduces  the  work  of  both  analysis  and 
computation  considerably. 

The  earliest  STC  is  Kalman's  [51]  self-optimizing  control,  which 
tried  to  adapt  the  controller  automatically  to  the  changes  of  the 
process  characteristics.  Peterka  [79]  extended  this  idea  to  the 
stochastic  system  with  constant  parameters.  It  is  Astrom  and 
Wittenmark  [10]  who  first  analyzed  the  properties  of  a  self-tuning 
regulator,  and  showed  that  optimal  control  can  still  be  obtained 
despite  the  approximation  introduced  by  the  certainty  equivalence 
principle.  Since  then,  the  interest  in  STC  increased  rapidly. 
Further  analysis  of  convergence  and  stability  has  provided  much 
deeper  insight  into  the  STC  theory.  Various  algorithms  for  both 
identification  and  control  have  been  proposed  to  make  STC  more 
efficient  and  versatile.  And  many  successful  applications  of  STC  have 
been  reported  justifying  its  substitution  for  PID  controllers. 
Furthermore,  the  sampled  data  based  algorithms  of  STC  are 
particularly  suitable  for  today's  low  cost  microprocessor  control. 


Most  of  these  developments  can  be  reviewed  by  studying  the  following 
representative  work.:  Astrom  [6,7];  Astrom  and  Wittenmark  [12];  Clarke 
and  Gawthrop  [27,28];  Goodwin  and  Sin  [39];  Harris  and  Billings  [41]; 
Isermann  [45,46];  Unbehauen  [100];  Well6tead  [105]. 

3.2  Self-tuning  Techniques 

Assume  that  the  system  to  be  controlled  can  be  described  by  an 
ARMAX  model  [39] 


A(z  1)yt  «  B(*"X)ut_k  +  C(z-1)et  (3.1) 

where  k,  an  integer,  is  the  time  delay  of  the  control  system  and 

, .  -1  -1  "’'a 

A(z  )  =  1  f  al2  +  ...  +  an  z 


B(z  1}  =  b0  +  b  1Z_1  +  +  V2  b 


C(z  S  =  1  +  clZ_i  +  ...  +  cn  z 

c 

are  associated  with  the  dynamics  of  the  system,  y  ,  ut  and  et  denote 
the  system  output,  the  control  input,  and  the  white  noise  disturbance 
at  the  sampling  time  instant  t,  respectively. 

The  general  form  of  a  self-tuning  controller  can  be  expressed  as 

=  F(z  )y.  +  G(z  1 )u_ 


-1  * 

HU  >y_,. 


(3.2) 


m 


where  y  is  the  desired  output  of  the  system.  The  polynomials  F,  G 


and  H  are  determined  by  the  polynomials  A,  B  and  C  in  (3.1)  and  the 


specified  control  criteria.  Depending  on  the  methods  for  determining 


the  parameters  of  the  controller,  self-tuning  controllers  can  be 


divided  into  two  categories,  explicit  and  implicit  algorithms.  The 


explicit  or  indirect  algorithm  needs  to  estimate  the  polynomials  A,  B 


and  C  before  assessing  the  parameters  of  the  controller.  In  the 


implicit  or  direct  algorithm  the  model  (3.1)  can  be  reparameterized 


into  a  certain  form  so  that  the  parameters  of  the  controller  can  be 


estimated  directly. 


3.2.1  Parameter  Estimation 


Both  the  explicit  and  implicit  algorithms  need  a  recursive 


parameter  estimation  scheme.  The  scheme  can  be  one  of  the  following 


recursive  estimation  methods  [ 46 , 33 ) : 


recursive  least  squares  (RLS) 


recursive  extended  least  squares  (RELS) 


recursive  maximum  likelihood  (RML) 


recursive  instrumental  variables  (RIV) 


stochastic  approximation  (STA) 


The  RML  or  recursive  prediction  error  method  (RPEM)  is  the  most 


general  and  complex  algorithm  [59],  It  has  a  unified  form  as  shown 


below: 


®(t)  “  O(t-l)  +  L(t)e(t) 


(3.3a) 


L(t)  =  P(t)i|/(t) 


1  +  *T(t)P(t-l)*(t) 

P(t)  -  [I  -  L(t)/(t)]P(t-l)  /  a(t) 

e(t)  -  y(t)  -  ?(t) 

-  y(t)  -  9  (t-l)4>(c) 


-  [5  (t),...,an  (t),U1(t),...,Bn  (t), 
a  b 


e  (t),...,e  (t ) ] 

1  n 
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+  2  (t-l)z  )  u(t) 
n 


I  Ci(t) 
E(t) 


(3.4e) 


where  matrix  I  is  a  unit  matrix  and  o(t)  is  the  forgetting  factor  for 

discounting  the  old  data.  The  a  posteriori  prediction  error  e(t)  in 

(3.4b  and  c)  can  improve  the  convergence  rate  of  the  estimation 

algorithm.  Most  of  the  time,  however,  it  is  approximated  by  e(t)  for 

simplicity.  Further  simplification  can  still  be  made.  If  <Kt)  is 

approximated  by  4>(t),  the  algorithm  reduces  to  RELS.  If  n  is  equal 

c 

to  zero,  i.e.,  C(z  )  ■  1,  the  algorithm  reduces  to  RLS.  RIV  and  STA 
can  be  obtained  by  modifying  or  simplifying  RLS;  however,  they  are 
seldom  used  in  STC. 


3.2.2  Control  Criteria 


In  general,  there  are  two  types  of  control  criteria,  classical 
and  optimal  [109].  From  (3.1)  and  (3.2),  the  closed  loop  system  is 
of  the  form 


(AG+BF)y  -  BHy  ,  +  CGe  , 
t+k  1  t+k  t+k 


(3.5) 


With  the  classical  type  ,  the  parameters  of  the  controller,  F  and  G, 
are  determined  such  that 


AA  a 

AG  +  BF  -  CT 


(3.6) 


where  A,  B,  and  C  are  the  estimated  polynomials  of  A,  B,  and  C.  T  Is 

the  prespecified  characteristic  polynomial  of  the  closed  loop  system. 

* 

For  a  regulator  system,  y  ■  0,  T  is  chosen  to  have  the  properties  of 
good  noise  rejection,  a  pole  placement  control  [106,107].  For  a 
servo  system,  et  =0,  T  and  H  are  designed  to  assure  good  tracking 
capability.  It  belongs  to  the  pole-zero  placement  control  [11].  The 
extended  STC  of  Wellstead  and  Sanoff's  [108]  chooses  T  and  H  for  both 
output  tracking  and  noise  rejection. 

With  the  optimal  type,  the  controller  of  (3.2)  is  taken  to 
minimize  the  cost  function 

I  *  EU2(t+k)|t] 

<Kt+k)  -  Pyt+k  +  Qut  -  Ryt+k  (3.7) 

where  P,  Q  and  R  are  polynomials  in  operator  z  Consequently, 

A  ^  -A 

F ,  G  and  H  can  be  determined  from 

~  ~~  -k- 

PC  -  AE  +  z  F 

G  «  BE  +  QC 

H  -  CR  (3.8) 

The  alternative  form  of  (3.5)  is  thus 
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This  is  called  generalized  minimum  variance  control  [23,25,26,38]* 

* 

When  P  ■  1,  y  ■  0,  and  Q  ■  0,  the  controller  reduces  to  the  famous 
self-tuning  regulator,  a  minimum  variance  controller  l 10].  Normally, 
the  cost  function  is  specified  with 

P  -  R  =  1,  Q  ■  A  >  0 
for  weighted  minimum  variance  control  [39]. 

The  cost  function  of  the  optimal  type  can  be  extended  to  state 
space  variables  for  more  sophisticated  control  156];  however,  it  is 
beyond  our  scope  of  interest.  The  generalized  STC  proposed  by 
Allidina  et  al.  [3]  is  a  unification  for  the  classical  type  and  the 
optimal  type. 

3.3  Implementation  Aspects 

In  the  practical  implementation  of  STC  there  are  several  factors 
which  should  be  considered.  They  can  be  classified  into 
identification,  control,  and  general  issues. 

3.3.1  Identification  Issues 

In  Identification  or  parameter  estimation,  we  consider  the 
following  factors: 

i.  Orders  of  system,  n  ,  n,  ,  and  n  :  They  should  be  known  a 

a  d  c 

priori.  Normally,  by  knowing  the  upper  bound,  we  can  let 

n  -  n.  *  n  -  n.  For  implicit  6elf-tuning  algorithms,  the 
a  b  c 

information  on  time  delay  k  is  also  needed.  These  integers 
determine  the  total  number  of  parameters  to  be  estimated. 


ii.  Bias  of  estimation:  Explicit  algorithms  require  identification 
schemes  with  less  bias  error  in  estimated  parameters,  such  as 
RML  and  RELS.  These  schemes,  however,  are  complex 
computationally  and  are  susceptible  to  numerical  instabilities. 
Implicit  algorithms,  on  the  other  hand,  may  use  the  simple 
scheme  as  RLS  and  still  have  satisfactory  control  performance. 

iii.  Computation  accuracy  and  speed:  In  microprocessor  based  control 

systems  the  accuracy  of  data  words  is  much  more  limited  than 
that  in  general  purpose  computers.  The  accumulated 
computational  errors  may  cause  numerical  difficulties.  To 

overcome  this  problem,  the  square  roots  method  [92]  and  the  UD 
method  [14,96]  of  RLS  can  be  used  to  convert  the  computation 
from  single  precision  to  double  precision  without  slowing  down 
the  computational  speed.  If  the  computational  speed  is 
critical,  a  fast  algorithm  is  recommended  [71,57]. 

iv.  Forgetting  factors:  A  forgetting  factor  is  very  important  in 
STC.  It  not  only  determines  the  tracking  speed  of 
identification  for  systems  with  time  varying  parameters,  but 
also  may  affect  the  convergence  of  the  STC  output.  Normally,  it 

is  in  the  range  of  .95  to  1.0.  A  low  forgetting  factor  is  good 
for  initial  coarse  tunings,  fast  varying  parameters,  and  abrupt 
changes  in  the  desired  output.  However,  it  may  destabilize  the 
control  system  if  the  value  is  too  low,  and  may  cause  estimator 
wind-up  in  servo  control  [5,12].  A  time  varying  forgetting 
factor  [36]  can  be  used  to  avoid  those  problems. 


v.  Initial  data:  This  is  not  a  crucial  issue  in  STC.  The  general 
rule  is  to  make  an  initial  guess  of  the  parameters  as  close  to 
ti.e  true  values  as  possible.  If  the  confidence  of  the  guess  is 
low,  the  determinant  of  the  covariance  matrix  should  be  large. 

vi.  Identif iability:  Astrom  and  Wittenmark  [10]  showed  that  in  an 
implicit  self-tuning  regulator  it  is  possible  to  have  the 
problem  of  identif iability  in  estimation.  To  avoid  any 
numerical  difficulties,  the  simplest  way  is  to  give  a  certain 
parameter  a  constant  value. 

3.3.2  Control  Issues 

The  selection  of  control  criterion  depends  on  the  following 

factors: 

i.  Performance  specifications:  This  includes  steady  state  error, 
noise  rejection,  convergence,  stability,  and  robustness.  These 
factors  usually  conflict  with  each  other.  Therefore,  the 
selected  criterion  must  be  a  compromise  of  these  factors.  At 
this  moment  we  are  unable  to  evaluate  or  justify  any  control 
criteria  for  the  lack  of  theoretical  analysis.  Some  general 
guidelines  can  still  be  found  [12]. 

ii.  Nonminimura  phase  problem:  The  nonminimum  phase  problem  refers 
particularly  to  some  of  the  zeros  of  B  in  (3.1)  falling  outside 
the  unit  circle  Minimum  variance  control  is  very  likely  to  be 
unstable  for  this  kind  of  systems.  Pole  placement  or  weighted 
minimum  variance  control  can  overcome  this  problem  [24]. 


iii.  Constraints  on  control  input:  Another  deficiency  of  minimum 

variance  control  is  the  large  amplitude  or  energy  of  control 
inputs  frequently  required  during  control.  The  control  hardware 
or  actuator  can  be  easily  overloaded  or  saturated.  Like  the 
nonminimum  phase  problem,  pole  placement  and  weighted  minimum 
variance  control  are  possible  solutions.  Or,  the  constraints 
can  be  incorporated  into  the  cost  function  of  optimal  control 
[62,61]. 

iv.  Deterministic  disturbances:  In  this  text,  deterministic 

disturbances  refer  to  the  slowly  varying  disturbances  and  the 
periodic  disturbances  with  unknown  frequencies  and  amplitudes. 
This  type  of  disturbances  is  important  in  practical  control 
systems;  however,  they  are  rarely  discussed  previously  in  STC. 
The  analysis  and  simulation  of  Goodwin  and  Sin  [39]  reveal  that 
the  deterministic  disturbances  can  be  well  compensated  for  by 
STC. 

3.3.3  General  Issues 

i.  Simplicity  in  implementation:  STC  should  always  be  designed  as 
simple  as  possible  for  practical  applications.  The  self-tuning 
algorithms  with  weighted  minimum  variance  criterion  and  RLS 
estimation  usually  give  more  satisfactory  results  than  other 
sophisticated  algorithms.  Also,  implicit  algorithms  are  simpler 
and  more  robust  than  explicit  ones  in  general  [12]. 

ii.  Sampling  rate:  Sampling  rate,  in  general,  should  be  high  enough 
to  cover  the  desired  bandwidth  of  the  closed  loop  system.  Its 


upper  bound  is  limited  by  some  considerations.  It  should  not 
exceed  the  response  speed  of  the  actuator.  It  should  be  low 
enough  so  that  the  computer  can  accomplish  all  the  digital 
signal  processing,  and  the  variation  of  system  delay  is 
insignificant  compared  to  the  sampling  time  interval.  Further 
considerations  can  be  seen  in  [12]. 
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CHAPTER  4 


THEORETICAL  DEVELOPMENT 

As  mentioned  in  Chapter  2,  the  geometrical  inaccuracies  of  the 
machined  part  result  from  the  relative  error  motion  between  the  tool 
tip  and  the  part  surface  during  machining  operations.  The  error 
motion  in  radial  direction  is  particularly  important  in  determining 
the  geometrical  errors.  There  are  numerous  sources  which  contribute 
to  such  kind  of  error  motion.  Three  major  types  of  error  sources 
considered  in  this  text  are:  stochastic,  deterministic,  and 
regenerative  types. 

The  stochastic  type  is  mainly  induced  by  the  play  in  the 
imperfect  spindle  bearings.  Fig.  2.10  is  a  typical  polar  plot  of 
spindle  error  motion  measured  from  the  idling  spindle.  The 
randomness  of  the  error  motion  suggests  statistical  models,  such  as 
autoregressive  model  [83]  and  moving  average  model,  for  describing 
this  error  motion  mathematically.  In  real  cutting,  this  error  motion 
is  complicated  further  by  the  chip  breakage,  nonstationary  built-up- 
edge,  hard  grains  in  work  material,  floor  noises,  and  change  of  the 
bearing  play  due  to  cutting  forces  and  thermal  drifts. 

The  deterministic  type  can  be  further  divided  into  linear  errors 
and  periodic  errors.  Linear  errors  refer  to  those  caused  by  steady 
or  slowly  varying  disturbances,  such  as  the  structural  deflection 


induced  by  the  mean  cutting  forces,  tool  wear,  and  thermal  expansion. 
Periodic  errors  mostly  come  from  the  unbalanced  rotating  parts,  such 
as  the  spindle,  grinding  wheel,  and  gear  train.  The  resulting  error 
motion  is  characterized  by  the  periodic  movement  with  frequencies 
having  integer  multiples  of  the  angular  frequency  of  the  spindle. 
The  fundamental  or  lowest  frequency  component  of  error  motion  results 
in  the  eccentricity  of  the  machined  part,  while  the  higher  harmonics 
produce  lobe-shape  errors  on  the  part  surface.  Periodic  errors  may 
also  result  from  directional  stiffness  variation  caused  by  the 
chuckhead  [81]. 


The  regenerative  type  is  the  regenerative  effect  discussed  in 
Chapter  2.  It  relates  to  the  interaction  between  the  cutting  force 
and  the  machine  tool  structure.  The  frequency  of  this  type  of  error 
motion  is  close  to  the  natural  frequency  of  the  tool-workpiece 
structure.  It  is  an  unstable  phenomenon  during  machining  operations. 
Whether  this  unacceptable  error  motion  will  occur  or  not  is 
determined  by  the  stiffness  and  the  damping  of  the  machining  system. 


Most  of  the  above  error  sources  are  dynamic.  They  are  difficult 
to  identify  because  the  error  sources  are  interrelated.  Moreover, 
these  error  sources  differ  with  part  designs  and  machining 
conditions.  Therefore,  they  can  be  best  compensated  for  by  a  GAC 
system.  The  GAC  system  proposed  in  the  next  section  assumes  that 
direct  measurement  of  the  part  geometry  is  used.  Thus,  the  effect  of 
the  above  major  error  sources,  except  the  thermal  effect,  on 
geometrical  accuracy  can  be  detected  and  corrected  by  the  controller. 
The  thermal  effect  is  the  dimensional  expansion  of  the  workpiece  due 


to  the  cutting  heat.  The  GAC  system  tends  to  generate  the 
geometrical  accuracies  at  the  elevated  workpiece  temperature  during 
machining.  Thus,  a  certain  temperature  measurement  and  calibration 
is  needed  for  the  true  dimensional  accuracy.  In  this  text,  however, 
this  effect  is  assumed  to  be  negligible  for  simplicity. 

4.1  Model  of  the  GAC  System 

The  physical  model  of  a  cylindrical  machining  process,  plunge 
cutting  in  particular,  is  shown  in  Fig.  4.1.  The  tool  cuts  into  the 
rotating  workpiece,  which  is  clamped  to  the  spindle.  The  forces 
generated  during  the  cutting  process  excite  the  tool-workpiece 
structure,  which  is  basical lv  a  spring-mass-damper  system. 
Meanwhile,  the  induced  structural  or  regenerative  type  error  motion 
together  with  the  other  two  types  of  error  motions  determine  the 
actual  tool  path  on  the  workpiece.  The  instantaneous  part  surface 
and  the  depth-of-cut  are  thus  formed. 

Based  on  the  above  investigation,  a  control  diagram  of  the  GAC 

system  for  plunge  cutting  is  constructed  as  shown  in  Fig.  4.2.  The 

profile  of  the  part  surface  at  the  tool  tip  is  denoted  by  ,  which 

is  the  reduction  of  radius  measured  from  a  base  circle  as  shown  in 

Fig.  4.1.  Due  to  the  regenerative  effect  in  plunge  cutting,  the 

instantaneous  depth-of-cut  is  Y  -Y  ,  where  pis  the  time  interval 

t  t-p 

for  the  spindle  to  rotate  in  one  revolution.  The  cutting  process  CP 
is  related  to  the  depth-of-cut  and  the  cutting  force.  Through  the 
structural  dynamics,  MT,  a  spring-mass-damper  system,  the  error 
motion  X  is  induced.  Thus,  surface  profile  Y  is  the  summation  of 
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Figure  4.1  Structural  Dynamics  and  Regenerative  Effect 
of  Plunge  Cutting 


Adaptive 

Controller 

Figure  4.2  Control  Diagram  of  the  GAC  System 


structural  error  motion  -X  ,  deterministic  error  motion  Z  , 

t  t 

stochastic  error  motion  W  ,  specified  infeed  of  the  tool  Y  (t) ,  and 

t  s 

the  control  action  u  from  the  adaptive  controller.  The  controller 

t-k 

G  constantly  measures  the  geometrical  error  y  ,  defined  by  Y  -Y  (t), 

t  t  s 

identifies  the  system  model,  and  gives  appropriate  control  action  to 

reduce  y  . 

t 

From  Appendix  1,  a  linear  model  of  the  GAC  system  is  derived  ant1 
expressed  by  a  stochastic  difference  equation  as  follows: 


A(  z  )  (y  -  y  +  <$  ) 
t  t-p  t 


"  B(z  1)(Vk  ■  VP  +  V  +  c(z  1)5t 


(4.1) 


-1  -1  -n 

A(z  )  =  1  +  3^  z  +  ...  +  anz 


-1  -1  -n 

B(z  )  =  b  +  b  z  ■*-...  +  b  z 
0  1  n 


(4.2) 


-1.  ,  -1  -n 

C(z  )  =  1  +  c  z  +  ...  +  c  z 
1  n 

where  z  *  is  the  backward  shift  operator,  that  is,  z  y  -  yt_j.  The 

inclusion  of  y  in  (4.1)  explains  the  regenerative  effect  in  plunge 
t-p 

cutting  process.  6  ,  equal  to  Y  (t)  -  Y  (t-p),  is  the  d  pth-of-cut 

t  s  s 

corresponding  to  the  specified  infeed.  The  coefficients  of  A(z  ) 


and  B(z  )  are  parameters  associated  with  the  characteristic 

-1 

properties  of  CP  and  MT.  They  are  not  known  in  general.  C( z  K 


characterizes  the  stochastic  error  source  in  machining,  where  C 

t 

2 

represents  a  white  noise  process  with  zero  mean  and  variance  o^.  The 
deterministic  error  source  is  embedded  in  (4.1)  by  the  common  zeros 
of  A,  B,  and  C  on  the  unit  circle.  In  the  following  control 
algorithms  the  order  of  the  model,  n,  and  time  delay,  k,  are  assumed 
to  be  known.  The  delay,  k,  which  accounts  for  the  computation  or 
measurement  delay,  is  far  smaller  than  the  regenerative  delay  p  in 
practice.  The  minimum  value  of  k  in  digital  control  system  is  equal 
to  one. 

4.2  STC  Algorithm 

By  solving  the  Diophantine  equation 

C(z_1)  *  ACz'SeCz'1)  +  z~kF(z~1)  (4.3) 


where 


-lx  -1  -k+1 

E(z  )  =  1  +  e  z  +  ...  +  e  z 
1  k-1 


-1  -1  -n+1 

F(z  )  -  fQ  +  fj*  +  ...  +  VlZ 


Equation  (4.1)  can  be  written  as 


C(z  1 )y  =  F(z  1 )(y  -  y  +  <5  )  + 

t  t-k  t-p-k  t-k 


G(z  S(u  -  y  +  6  )  + 
t-k  t-p  t 


C(z  J)(y  -  6  )  +  C(z  *)E(z  Sc 
t-p  t  t 


(4.4) 
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where 


G(z  * )  =  B(z  *  )E(z  S 


Thus, 


F  G 

y  =  — (y  -  y  +  6  )  +  — (u  -  y  +6  )  + 

yt+k  C  yt  yt-p  t  C  t  't-p+k  t+k 


(y  -6  )  +  E£ 

t-p+k  t+k  t+k 


yt+k|t  +  E^t+k 


(A. 5) 


-1 

where  the  operator  z  is  omitted  for  simplifying  the  representation. 

y  is  the  conditional  expectation  of  y  based  on  the 

t+k I t  t+k 


observations  y  ,  y 

t  t-1 


and  the  parameters  of  A,  B, 
and  C.  Thus,  for  the  minimum  variance  control  criterion 


u  ,  u 
t  t-1 


Min  E [ y 


2 

t  +  k  1 


Ei  .  J :  an  exp.  tat  ion  t unction 


(  A  .  b  ) 


it  gives  the  following  control  law 


V  t  +  k|  t 


U 


(A. 7) 


or 


Cy 


t  +  k  t 


f(y. 


y  +  6  )  +  G(u 
t-p  t  t 


y  + 

t-p+k 


6 

t+k 


)  + 


C(y  -  6  )  -  0 

t-p+k  t+k 


(A. 8) 


with  variance 


<o 


(4.9) 


Var  [ y t ]  ■  [i  +  +  ...  +  e^Jo 

In  self-tuning  control,  since  all  the  parameters  are  unknown,  (4.8) 
should  be  written  as 


9  .  +29  i  +...+29  i 

r t+k 1 1  r  t+k— 1 1 1— 1  n^t+k-n|t-n 


F(y  -  y  +  6  )  +  G(u  -  y  +6  )  + 

t  t-p  t  t  t-p+k  t+k 


C(y  -  6  ) 

t-p+k  t+k 


*  0  (4.10) 

where  all  the  estimated  quantities  are  represented  by  placing  the 

symbol,  %  on  the  top  of  the  notations.  Therefore,  based  on  the 

certainty  equivalence  principle,  the  control  law  corresponding  to 

(4.7)  is  to  find  u  such  that 
t 


?t+k|  t 


y  +  6  )  +  G(u 
t-p  t  t 


-  y  +6  ) 

t-p+k  t+k 


+  C(y 


t-p+k 


*  0  (4.11) 

(4.11)  is  derived  from  (4.10)  by  dropping  the  terms,  9  *  iMl, 

t+k-i | t-i 

A 

2,...,  n.  Consequently,  the  model  for  estimating  the  parameters,  f, 

A 

,  and  C,  is 
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y  *  F(y  -  y  -  6  )  +  G(u  -  y  +  6  ) 

t  t  t-p+k  t  k  t-k  t-p  t 

C(y  -  6  )  +  e  (4.12) 

t-p  t  t 

where  e  corresponds  to  EC  •  If  the  system  is  optimally  controlled, 
the  control  output  y  should  satisfy  (4.11). 

The  STC  algorithm,  represented  by  (4.11)  and  (4.12),  belongs  to 
the  implicit  type.  All  the  parameters  in  (4.12)  can  be  updated  by 
the  recursive  least  squares  (RLS)  scheme  as  shown  in  Appendix  2.  The 
recursive  scheme  can  be  further  subjected  to  UD  factorization  for 
improving  the  accuracy  of  computation  [27,14].  A  forgetting  factor 
is  used  for  tracking  the  specified  input  and  the  time-varying 
parameters.  Also,  it  is  helpful  to  enhance  the  convergence  rate  of 
the  output  at  the  outset  of  control. 

4.3  Stability,  Convergence,  and  Robustness 

In  optimal  control,  it  can  be  shown  that  the  control  law  (4.7) 
gives  a  closed  loop  system  satisfying 

Cyt  -CE£t 

CBu  -  C(B  -  A) (y  -  6  )  -  CF£  (4.13) 

t  t-p+k  t+k  t 

The  characteristics  of  the  closed  loop  system,  thus,  depend  on  the 
polynomials,  B  and  C.  That  is,  for  a  stable  minimum  variance 
control,  the  zeros  of  B( z  *)  and  C(z  S  must  be  within  the  unit 


circle 


The  zeros  of  A,  however,  plays  no  important  role  in  the 


minimum  variance  control  system.  Therefore,  even  though  the  open 
loop  system  is  unstab) e,  this  closed  loop  system  may  remain  stable 
and  optimally  controlled.  Conceivably,  the  self-tuning  control 
system  will  inherit  the  properties  of  the  optimal  control  system. 
This  can  be  seen  from  the  following  stability  and  convergence 
analysis  conducted  by  Goodwin  and  Sin  [39,88]. 

The  existence  of  C  in  (4.12)  is  due  to  the  specified  infeed  and 
the  regenerative  effect.  Hence,  it  is  a  servo  control  problem,  by 
introducing  the  following  substitutions 

y*(t)  -  Ys(t) 

it 

y(t)  -  yt  +  y  (t) 

* 

u(t-k)  -  uc_k  +  y  (t) 

A'  -  A  +  (B  -  A)z'P 

w(t)  -  (■  (t) 

the  model,  (4.1),  can  be  converted  into 

A'y(t)  -  Bu(t-k)  +  Oo(t) 
and  the  control  law,  (4.10),  into 

4>(t)  6(t)  *  y  (t+k) 
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where  <j>(t)  and  9(t)  are  defined  in  [39].  The  term  u>(t)  is  assumed  to 
be  a  white  noise  process,  which  satisfies  the  following  conditions: 

(i)  E[u>(t)|t-1]  «  0 
(ii)  E[to2(t)  1 1— 1  ]  «  a2 
1  N  2 

(iii)  lim  Sup  —  E  ui(t)  <  00 
t«l 

Also,  the  following  assumptions  are  made  about  the  system: 

1 .  k  is  known . 

2.  The  upper  bound  for  model  order  is  known. 

3.  B(z  *)  and  C(z  1)  have  all  zeros  inside  the  unit  circle. 

4.  £Yz7  _  \  is  strictly  passive. 

By  using  a  modified  RLS  [88] ,  the  self-tuning  controller  can  be  shown 
to  ensure  with  probability  one  that 


N 


lim  Sup  i  l  y(t)2  <  <*> 
Nx»  wt-l 


(4.14) 


N 


lim  Sup  t-t  1  u ( t ) 2  <  » 
t-1 


(4.13) 


lim  Sup  4  ^  E[(y(t)  -  y  ( t ) ) 2 1 1- 1  ] 


k-1 

2  „  2 
o  E  e. 


(4.16) 


Properties  (4.14)  and  (4.15)  yield  the  stability  of  the  closed 

loop  system.  The  energy  of  both  the  system  output  and  the  control 

input  are  bounded.  Condition  (4.16)  stands  for  the  output 

* 

convergence.  That  is,  the  desired  output  y  (t)  is  asymptotically  and 

optimally  tracked.  The  passivity  or  positive  real  assumption  on 
1  1 

~  -  ~  a  stronger  assumption  than  C  having  asymptotically  stable  zero 

v  4- 

[58],  is  important  to  the  output  convergence. 

From  Appendix  1,  A',  B,  and  C  may  have  common  zeros  on  the  unit 
circle  because  of  the  deterministic  disturbance.  When  the  sampling 
rate  is  sufficiently  high  B  may  have  unstable  zeros  due  to  the  excess 
number  of  poles  [8].  As  for  the  former  case,  the  common  zeros 
characterize  the  uncontrollable  part  of  the  system.  In  optimal 
control,  Goodwin  and  Sin  [39]  showed  that  if  C  is  split  into 

C  -  C  +  C 
s  r 

where  C  has  only  asymptotically  stable  zeros  and  C  is  the  remainder 
s  r 

which  can  be  as  small  as  desired  by  the  choice  of  C  ,  the  use  of  C 

s  s 

instead  of  C  in  the  Diophantine  equation  (4.3)  will  result  in  a 
suboptimal  but  stable  control.  By  the  same  token,  in  self-tuning 
control,  if  the  zeros  of  $  is  kept  within  the  unit  circle  by 
projection  method  [39],  the  closed  loop  system  will  be  nearly  optimal 
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For  Che  later  case,  u  in  (4.13)  will  be  unbounded  when  B  has 

t 

unstable  zeros.  If  lowering  the  sampling  rate  is  not  permissible, 
the  weighted  minimum  variance  control  as  shown  in  Appendix  2  can  be 
used  to  get  a  suboptimal  but  stable  control. 

The  above  analysis  is  based  on  the  assumption  that  the  model 
order  n  is  sufficiently  large  to  cover  all  the  natural  frequencies  of 
the  system.  In  practice,  a  lower  order  model,  however,  is  desirable 
to  simplify  the  algorithm  and  minimize  the  computation  time.  This 
prompts  the  robustness  issue  of  self-tuning  control.  Again,  the 
weighted  minimum  variance  control  can  be  used  for  ensuring  the 
robustness  of  the  low  order  controller  [12J.  Furthermore,  by 
properly  varying  the  weighting  factor  on-line,  the  control  input  can 
be  assured  to  fall  within  the  constraints. 

4.4  Multiprobe  Measurement  in  GAC 


To  separate  the  spindle  error  motion  from  the  measurement  of  the 
surface  profile,  Whitehouse  [111]  devised  the  multiprobe  method, 
which  uses  three  probes  of  different  sensitivities  without  requiring 
any  masters.  The  probes  are  spaced  around  the  workpiece  as  shown  in 
Fig.  4.3.  The  centerlines  of  these  probes  intersect  at  a  fixed  point 
0.  Due  to  the  spindle  error  motion,  the  axis  of  rotation  0'  moves 
around  the  fixed  center  randomly.  The  amplitude  of  this  movement,  e, 

is  assumed  to  be  far  less  than  the  nominal  radius  of  the  part, 

-3 

normally  below  10  R  [110]. 


6o 


y.(t)  =  R  -  r.(t) 

V.(c)  =  R  -  d  .(c)  ,  i  =  1,  2,  3.  (4.17) 

where  r  (c)  and  d  (t)  are  illustrated  in  Fig.  4.3.  Hence, 
i  i 

v.(t)  =  y .  (t)  +  ( r  ( t )  -  d  (t)),  i  =  1,  2,  3.  (4.18) 

ii  l  i 

It  can  be  seen  from  Fig.  4.3  and  (4.18)  that  the  spindle  error  motion 

involves  two  unknown  quantities:  amplitude,  e,  and  angle,  <J>.  At 

least  three  probes  are  thus  needed  to  remove  the  spindle  error  motion 

from  the  probe  measurements;  however,  the  geometrical  error,  y  (t), 

i 

can  never  be  obtained  for  the  number  of  unknown  quantities  always 
exceeds  that  of  known  quantities  in  (4.18).  Therefore,  only  three 
probes  are  sufficient  for  solving  the  separation  problem.  .Additional 
probes  can  be  regarded  as  redundant. 

With  the  three  probes  as  shown  in  Fig.  4.3,  it  can  be  proved 

t  hat 

d  =  cv  (t)  +  av  (t)  +  bv  (t) 
t  1  2  3 

=  cy  (t)  +  ay  (t)  +  by  (t)  (4.19) 

1  2  3 

with 

a  -  sin  3  /  sin(  a  +  3) 

b  =  sin  a  /  sin(a  +  3)  (4.2'J) 


c 


-1 
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and 

a  +  6  *  it,  2n 
a  ,  d  >  0 

where  d  is  the  combined  measurement  of  probes  1,  2,  and  3  with  the 
t 

scaling  factors  -1,  a,  and  b,  respectively.  Thus,  (4.19)  shows  that 

the  spindle  error  motion  hc.s  been  separated  from  the  measurement  of 

the  surface  geometry.  During  the  plunge  cutting  process,  R  in  (4.17) 

should  be  replaced  by  R  .  Nevertheless,  (4.19)  and  (4.20)  hold  in 

this  situation.  Instead  of  applying  Fourier  analysis  to  reconstruct 

the  surface  geometry  as  Whitehouse  did,  the  combined  signal  d  is 

t 

used  in  the  GAC  system  by  the  following  manipulation. 

According  to  the  angular  spacing  of  these  probes  and  the 

direction  of  rotation,  the  geometrical  errors  y^(t),  y?(t),  and  y^(t) 

can  be  expressed  by  sampled  variables,  y  ,  y  ,  and  v 

t-f  t-'g  '  r 

respectively.  Thus,  (4.19)  can  be  written  as 

d  =  by  -y  +  ay  ,  t  *  1,  2,  3,... 
t  t  t-f  C-g 

f  =*  p  3/ 2  TT 

g  •  p( a  +  6)/ 2n 

where  p  is  the  number  of  sampled  points  per  revolution.  Define 

-1  — f  -g 

P(z)*b-z  +  az  (4.21) 

d 


dt  =  pd(2  )yt 


(4.22) 


Before  proceeding  to  the  design  of  a  GAC  system,  the  characteristics 


of  the  polynomial,  P  (z  ),  should  be  analyzed. 

d 


4.4.1  Properties  of  Multiprobe  Measurement 


The  characteristic  roots  of  P  (z  )  can  be  found  by  solving 

d 


Pd(z  1 )  *  0 


g  g-f 

z  sin  a  -  z  sin(  a  +  3)  +  sin  3=0  (4.23; 


j  2  irra/  p 


,  j  -  v/-T 


where  m  is  a  real  number  and 


0  <  m  <  p 


then  (4.23)  can  written  as 


r  [cos  m( a  +  3)  +  j  sin  m(  a  +  3)] 


CT—  f 

-  r  (cos  ma  +  j  sin  ma)sin(a  +  3)  +  sin  3=0 


v  ■  .•  •* 


.**  .**  .*• 


•  »  •  •  „*  .*  ,*  ,=  mm  ,  „ 

*s- Y  sJCj-m i.  *u  ^  vC'vV WV*" Jv' 


And,  in  this  case,  (4.24)  gives 


(±)r8sin  a  -  (±)r8  sin(a  +  0)  +  sin  0=0  (4.28) 

*  * 

Therefore,  if  there  exists  a  characteristic  pair  (r  ,m  ),  it  roust 

★ 

either  satisfy  (4.26)  and  (4.27),  or  (4.28).  The  magnitude  of  r  can 
be  any  positive  finite  value.  One  special  characteristic  pair  is 


*  * 
r  =  1 ,  m  =  1 


for  all  admissible  a  and  0.  The  solution  means  that  if  y  is  a 

t 

sinusoidal  function  with  frequency  — 


then 


=  sin  2itt/p 


P  (z  1  )y  =0 
d  t 

In  measurement,  this  means  that  the  multiprobe  method  can  not  detect 
the  eccentricity  error  of  the  part,  that  is,  the  sinusoidal  form 
error  with  once-per-revolution  frequency.  This  is  the  major 
constraint  of  multiprobe  method  [111],  If  this  error  is  caused  by 
the  unbalanced  motion  of  a  certain  machine  member,  a  balancing 
operation  should  be  performed  before  machining.  The  higher 

*  it 

harmonics,  corresponding  to  r  =1  and  m  =  2,3,4,...,  may  also  be 
undetectable  to  the  measurement  setup  if  a  and  0  are  not  properly 
chosen.  To  avoid  this,  the  following  analysis  is  needed. 


occurs 


From  (4.26)  and  (4.27),  the  suppression  of  higher  harmonics 
* 

when  m  *  m  and 


m  m  /p  ■  k  (1  ±  m  /p)  -  1,  2,  3, 
a  a  a  ’ 


m  m0/p  *  k3  (1  1  me/p)  "  l»  2»  3‘ 


or 


*  k  P  kQP 

*  a  3 

m  -  m  ±  1 ,  m  »  -  -  - 

0  0  m  m„ 

a  B 


(4.29) 


From  (4.29) 


m  k  ■  m  k  -  k[lcm(m  ,m  )) 
a  B  B  a  a  B 


where  "lcra"  is  lowest  common  multiplier.  And 


kp[lcm(m  ,ra  )] 
_ a  6 

m  m  n 
or  p 


_ kp. . . . 

gcd(ma,mg) 


where  "gcd "  is  greatest  common  divisor  and  k  can  be  any  integers 
which  make  m^  an  integer.  Hence, 


gcd(ma,me) 


or  more  precisely, 


* 

m 


mod[' 


kp 


gcd(ma,me) 


*  1,  Pi 


(4.30) 


From  (4.30),  if  m  is  relatively  prime  to 
a 

relatively  prime  to  p,  then 


m  ,  or  gcd(m  ,m  )  is 
8  a  e 


* 

m  ■  1,  p-1  (or  -1) 

that  is,  no  higher  harmonics  are  suppressed.  Sometimes,  it  is 

neither  easy  nor  necessary  to  achieve  this  property.  A  probe 

* 

configuration  with  sufficiently  large  m  is  acceptable.  This  is 
because  of  the  fact  that  higher  harmonics  usually  have  very  low 
amplitudes. 

The  other  desired  property  of  the  multiprobe  method,  which  will 

be  important  to  the  following  controller  design,  is  to  select  angles 

* 

a  and  6  such  that  the  number  of  characteristic  pairs  having  r  >  1  is 
minimal  or  none. 


It  seems  that  we  are  not  able  to  derive  any  explicit  function  in 
terras  of  a  and  6  for  r  from  (4.26),  (4.27),  and  (4.28).  However,  the 
following  rules  can  help  us  approach  the  desired  property. 

Rule  1:  It  can  be  shown  that  if  a  equals  to  6,  the  characteristic 
pair  will  be 


* 

[r 


1,  ( 


k£ 


m 


1)] 


k  *  0,  1,  2,...,m-l 
a 


■  •'  --A--  .--v-  .v  .vy.v-  -A  .. 


Therefore,  if  m  is  relatively  prime  to  p,  there  are  no  zeros  of 
a 

(4.23)  outside  the  unit  circle,  and  there  is  no  suppression  of 
higher  harmonics  of  y^. 

Rule  2:  From  (4.23),  if  |sinB/sina|  <  1,  the  chance  to  get  all 
the  zeros  to  lie  within  unit  circle  is  high,  but  it  is  still  not 
guaranteed. 

After  this  preliminary  selection  process,  all  zeros  of  (4.23)  can  be 
checked  by  numerical  solution.  Reselections  of  a  and  8  may  be 
needed. 


4.4.2  STC  With  Multiprobe  Measurement 


Multiplying  (4.1)  by  p  ,  the  model  of  the  GAC  system  becomes 

d 


-1 


A(d  -  d  +  6/)  -  B(u'  -  d^  +  6/)  +  C(z  )E  '  (4.31) 

t  t-p  t  t-k  t-p  t  t 


where  d  is  the  measurement  from  the  multiprobe  method  and 
t 


-1 


u'  ,  •  P,(z  )u  . 
t-k  d  t-k 


(4.32) 


-1 


V  "  Pd<2  >6t 


-1 


5  '  -  PXz  K 

t  d  t 


For  the  minimum  variance  control  criterion. 


Min  E(dt+k|t] 


the  self-tuning  control  algorithm  is  similar  to  (4.11)  and  (4.12). 


The  model  for  parameter  estimation  Is 


d  «  F(d  -  d  +  6'  )  + 

t  t-k  t— p— k  t— k 


G(u  '  -  d  +6'  ,  )  +  C(d  -  6  ')  +  e  '  (4.33) 

t  t-p+k  t+k  t-p  t  t 


and  the  control  law  Is 


F(d  -  d  +  6  ')  +  G(u  '  -  d  ,  +  <$'  ) 

t  t-p  t  '  v  t  t-p+k  t+k7 


+  C(d  -  6'  )  ■  0 

t-p+k  t+k 


(4.34) 


Since  P  has  at  least  one  zero  on  the  unit  circle,  and  possibly  has 
d 


unstable  zeros,  the  control  input  uf  can  not  be  calculated  directly 


by  using  (4.32).  Any  errors  in  u '  and  will  not  be 


-1 


damped  out  by  .  The  solution  for  this  problem  is  to  compute 


("t'  +  “t-f  -  “t-g*  '  b 


and 


u 

t 


PU 


t 


where  p,  called  detuning  factor,  is  a  scalar  and 


(4.33) 


0  <  p  <  1 

The  introduction  of  u^  and  p  is  similar  to  that  used  in  [26] 


criterion  of  p  is  that  the  more  unstable  the  zeros  of  P  are,  the 

d 

smaller  p  is.  Normally,  a  smaller  detuning  factor  results  in  less 
optimal  output. 


The  previously  mentioned  properties  of  stability,  convergence, 
and  robustness  can  be  readily  applied  here.  Conceivably,  however, 
the  performance  would  be  slightly  affected,  since  the  objective 
function  is  d^  instead  of  y^,  and  the  detuning  factor  is  less  than 
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CHAPTER  5 


SIMULATION  OF  THE  GAC  SYSTEM 


Simulation  is  an  important  tool  for  evaluating  the  effectiveness 
of  the  self-tuning  control  algorithm.  The  stability  and  convergence 
properties  described  in  the  last  chapter  are  based  on  certain 
idealized  assumptions.  They  provide,  at  best,  qualitative  guidelines 
for  assuring  or  improving  the  algorithms.  Still,  there  are  some 
questions  in  real  applications  that  are  very  difficult  to  answer 
through  theoretical  analysis.  Most  of  them  are  related  to  the 
robustness  of  the  algorithms,  such  as:  nonlinearities  in  the  physical 
process  and  controller,  insufficient  number  of  order  in  the  system 
model,  finite  precision  in  numerical  computation,  and  so  forth. 
Simulation  can,  however,  easily  and  quantitatively  show  the  best 
performance  that  the  control  algorithm  can  achieve  under  these 
practical  considerations.  And  simulation  can  determine  the 
appropriate  numerical  range  for  the  ad  hoc  factors,  such  as:  the 
forgetting  factor  in  parameter  estimation,  the  weighting  factor  in 
weighted  minimum  variance  control,  and  detuning  factor  in  GAC  system 
with  multiprobe  measurement.  Also,  another  important  advantage  of 
simulation  is  that  It  can  provide  experience  in  significantly 
lowering  the  possibility  of  failures  or  damagf  .n  setting  up  a  real 


system 


There  are  two  approaches  for  simulating  a  physical  process: 
digital  and  analog.  The  analog  approach  is  to  simulate  the  process 
by  using  electric  or  electronic  wirings.  The  analog  one  is 
convincing,  while  the  digital  approach,  by  using  a  digital  computer, 
is  superior  as  it  is  versatile  enough  to  simulate  various  complex 
situations  in  the  real  system  such  as,  time  delays.  Therefore,  this 
chapter  is  devoted  to  the  simulation  on  a  general  purpose  digital 
computer,  a  CDC  6000  system.  The  simulation  with  analog  wirings  and 
microprocessor  controller  will  be  discussed  in  the  next  chapter. 

A  proper  simulation  model  is  not  merely  the  one  which  contains 
as  much  complexity  as  the  process  may  have,  but  it  must  also  be  one 
whose  dynamic  behavior  is  in  coincidence  with  the  practical 
observations.  For  machining  processes  it  is  well  known  that  chatter 
or  instability  will  occur  whenever  the  machining  system  has 
insufficient  stiffness  and  damping.  This  important  phenomenon  will 
be  shown  in  the  following  simulation. 

5.1  Simulation  Model  and  Conditions 

A  finishing  process  of  plunge  grinding  or  cutting,  as  shown  in 
Fig.  4.1  and  Fig.  5.1,  is  simulated.  The  specified  infeed,  Yg(t), 
increases  at  a  constant  rate  and  then  stays  at  a  desired  part 
dimension  for  two  or  more  revolutions.  In  this  text  we  assume  that 
this  feed  motion  is  provided  by  a  conventional  hydraulic  or  electric 
drive.  The  control  action,  which  is  usually  small  and  of  high 
frequency,  can  be  generated  by  a  separate  precision  actuator  embedded 
in  the  tool  holder  or  in  the  spindle  housing.  The  piezeoelectric 
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translator  154,21,19)  is  an  excellent  device  for  implementing  the 
precision  actuator.  The  separation  of  the  control  motion  and  the 
feed  motion  is  advantageous  in  releasing  the  complicate  coordination 
between  the  basic  NC  functions  and  the  feedback  control.  Also,  the 
specially  designed  tool  holder,  with  the  kit  consisting  of  a  sensor 
and  a  microcontroller,  can  be  easily  reinstalled  in  any  other  similar 
machine  tools.  This  will  simplify  the  manufacturing  of  the  machine 
tool  and  the  maintenance  of  the  precision  actuator. 

To  simulate  the  machining  dynamics  on  the  digital  computer,  the 
discrete  transfer  function  of  the  machine  tool  structure  given  in 
Appendix  1  is  used.  A  nonlinear  process,  Yt-max{Yt',  as  shown 

in  Fig.  5.2,  is  added  to  account  for  the  loss  of  contact  which  may 
occur  in  finish  cutting.  Nevertheless,  the  derived  self-tuning 
control  algorithm,  based  on  the  linear  model,  will  be  used  for 
testing  its  robustness.  Unless  stated  otherwise,  the  simulation  uses 
the  following  data: 

damping  ratio  of  MT  (S)  »  .05 

natural  frequency  of  MT  (uj^)  -  9.55  hz  -  60  rad/sec 
speed  of  work  spindle  (ft)  ■  100  rpm 
sampling  period  (T)  -  10  msec 
regenerative  time  delay  (p)  »  1/  (ftT)  «  60 
total  number  of  samplings  (N)  »  900  (i.e.  15  revolutions) 
system  time  delay  (k)  -  1 
feedrate  *  .24  unit/rev  -  .004  unit/sample 
The  dimensionless  data  of  feedrate  is  chosen  only  for  numerical 
convenience.  Thus,  appropriate  sensitivities  of  the  transducers  and 


the  power  amplifiers  are  assumed  for  the  simulated  geometrical  errors 

and  control  inputs.  The  selection  of  natural  frequency,  (d^,  end 

sampling  period,  T,  is  done  according  to  the  rule  of  thumb,  to^T  in 

the  range  of  0.25-1,  or  6-25  samples  in  each  period  of  the  resonant 

response  of  the  closed  loop  system  [113].  It  can  be  easily  shown 

that  for  different  the  simulation  results  will  be  identical,  if 

u>  T  remains  the  same.  Therefore,  the  simulation  results  in  this 
n 

chapter  can  be  applied  to  systems  having  higher  natural  frequencies. 
The  system  delay,  k,  can  be  greater  than  one  if  the  positions  of 
measurement  and  control  are  not  at  the  same  spot.  The  delay  due  to 
computation  is  assumed  to  be  negligible. 

Another  important  parameter  not  mentioned  above  is  the  stiffness 
ratio,  (k^/k^),  discussed  in  Chapter  2.  The  larger  the  ratio  is,  the 
more  unstable  the  uncontrolled  machining  system  is.  Therefore,  this 
ratio  will  be  varied  to  test  the  control  performance  at  different 
stability  levels  of  the  machining  system. 

The  simulation  is  divided  into  two  major  parts.  The  first  part 
is  to  test  the  performance  of  the  control  strategy  by  assuming  that 
the  geometrical  error  y  can  be  measured  directly.  The  second  part 
is  to  include  the  reality  of  measurement  problem  and  to  test  the 
feasibility  of  the  GAC  system  with  multiprobe  measurement.  Table  5.1 
summarizes  the  simulation  conditions  for  both  parts. 

In  the  simulation  results,  the  geometrical  error  yt  is 
represented  by  three  conventional  terms:  dimensional  error,  which  is 
the  mean  value  of  y  ,  out-of-roundness ,  which  is  the  maximum  peak  to 
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peak  error  of  y  ,  and  roughness,  which  is  the  root  mean  square  of  y  . 

t  t 

Since  roughness  refers  to  the  high  frequency  error,  its  value  will 
not  be  counted  and  shown  in  deterministic  systems.  Notice  that  the 
definitions  of  these  terms  are  a  little  different  from  those  used  in 
practice.  However,  they  are  convenient  for  evaluating  the 
performance  of  the  simulated  GAC  system.  These  data  are  shown  with 
polar  plots  for  comparison.  Time  responses  are  plotted  to  check  the 
performance  of  both  controlled  and  uncontrolled  systems.  Negative 
value  of  the  responses  means  that  the  tool  is  away  from  the  workpiece 
with  respect  to  the  specified  tool  path.  The  plots  of  the  estimated 
parameters  are  used  for  examining  the  behavior  of  the  STC  algorithms. 

5.2  Simulation  With  Ideal  Measurement 

In  the  first  part  of  the  simulation,  we  begin  to  test  a  stable 
machining  system  without  stochastic  disturbance  in  order  to  see  how 
well  the  deterministic  disturbance  can  be  compensated  for  by  the 
self-tuning  controller.  Then,  some  stochastic  disturbance  is  added 
to  the  system  for  examining  the  noise  rejection  capability  of  the 
controller.  The  increase  in  stability  of  the  self-tuning  controller 
can  be  visualized  by  increasing  the  stiffness  ratio.  The  importance 
of  incorporating  the  physical  insight  into  the  mathematical  model  in 
Chapter  A  can  be  seen  by  comparing  the  control  performance  with  that 
of  the  controller  which  uses  the  general  self-tuning  algorithm  as 


described  in  Chapter  3 


5.2.1  Deterministic  Systems 


Fig.  5.3  depicts  a  simulation  run  with  the  periodic  disturbance 

having  frequency  2 ft  and  amplitude  0.004.  Stochastic  disturbance  is 

not  added  in  simulation.  It  is  a  stable  machining  system  with 

stiffness  ratio,  k  /k  =1.  The  polar  plots  show  that  dimensional 

c  m 

error  and  out-of-roundness  are  almost  eliminated  in  the  GAC  system. 
Moreover,  Fig.  5.3  shows  that  the  GAC  system  makes  perfect  tracking 
of  the  desired  tool  path  through  the  control  input  which  is  nearly 
equal  but  opposite  to  the  uncontrolled  response.  Because  of  this, 
the  time-consuming  dwell  cutting  or  spark-out  process  is  not  needed 
for  acquiring  satisfactory  accuracy.  Hence,  we  can  have  the  tool 
retract  from  the  workpiece  smoothly,  at  the  last  revolution  of  cut, 
by  letting  6  in  equations  (4.11)  and  (4.12)  linearly  decrease  as 
shown  in  Fig.  5.1b.  Ihe  simulated  control  input  also  demonstrates 
this  kind  of  arrangement. 

The  estimated  parameters  in  Fig.  5.3  quickly  converge  to  steady 
state  values,  although  the  adopted  model  order,  n,  is  two  instead  of 
four,  the  full  order  of  the  simulated  system.  By  substituting  those 
steady  state  estimates  into  (4.12),  the  behavior  of  the  controller 
can  be  explained.  Since  the  controlled  y  is  very  near  to  zero  in 
the  middle  of  cutting,  the  estimated  model  can  be  written  as 

-1 

(  1.37062  -  .53975z  )<5  + 

-1  -2 

(.86733  -  1.69675z  i-  ,86733z  )(ut_j  +  6t) 

-  (1  -  .5364z  ^  +  .40897z  )6^  =  0 
In  addition,  at  that  time,  6  is  a  constant 
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(1  -  1.9562912z  +  z  )u  j  -  .004376  (5.1) 

The  polynomial  on  the  left  hand  side  of  Che  above  equation  can  be 
expressed  as 


1  -  2(cos  2ir<uT)z  +  z 

which  represents  a  sinusoidal  wave  of  frequency  u.  Solving  for  u> 
gives 


<u  «  200  rpm  »  20 

which  is  the  frequency  of  the  periodic  disturbance  in  the  simulation. 
This  explains  the  usage  of  a  sinusoidal  control  input  in  Fig.  5.3  for 
countering  the  periodic  disturbance .  In  addition,  let 

u  -  u  +  u  (5.2) 

t  p  s 

where  u^  is  the  sinusoidal  component  and  ug  is  the  6tatic  component. 
Substituting  (5.2)  into  (5.1)  gives 

u  «  (.1)6  -  (k  /  k  )6 
s  cm 

The  above  solution  explains  why  the  self-tuning  controller  can 
automatically  compensate  for  the  static  deflection  during  the  cutting 


operation 


Because  of  Che  discrepancy  of  Che  model  order  and  the 
nonpersistent  excitation  of  the  output,  the  converged  parameters  are 
merely  a  set  of  estimates  which  satisfy  the  minimum  variance 
criterion  as  shown  above.  When  the  specified  depth-of-cut  begins  to 
change  at  the  end  of  the  cut,  it  can  be  seen  in  Fig.  5.3  that  these 
estimates  also  change  to  a  new  set  of  values,  accordingly,  to 
maintain  the  optimality. 

* 

A  time  varying  forgetting  factor  is  used  in  the  above  simulation 
to  improve  the  convergence  rate  of  the  estimates  and  the  output ;  it 
is  of  the  form  shown  below: 


where  t  is  the  sampling  time  instant  at  t  •  1,  2,  3,  ...  ,  and 

-  .995,  f2  -  .5 

This  forgetting  factor  is  low  at  the  start-up  of  tuning  so  that  the 
inaccuracy  of  the  initial  guesses  can  be  forgotten  as  fast  as 
possible.  The  forgetting  factor  then  increases  to  f^  asymptotically, 
f^  is  slightly  smaller  than  one  to  permit  the  self-tuning  controller 
to  adjust  itself  for  the  environmental  changes  or  modeling  errors. 
Fig.  5.4  shows  a  similar  simulation  run  with  constant  forgetting 
factor 


X  -  .995 


It  shows  that  the  controlled  response  and  the  estimated  parameters 
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converge  far  slower  than  they  do  in  the  simulation  run  shown  in  Fig. 
5.3.  A  rule  of  thumb,  gained  from  simulation  experience,  for 
selecting  f^  and  f^  is 


.99  <  fA  <  1.0 


>1  <  f2  <  1.0 


Faster  convergence  rate  means  that  the  GAC  system  can  take  larger 
depth-of-cut ,  shorten  the  production  time,  and  still  have  good 
accuracy. 


As  it  was  mentioned,  the  periodic  disturbance  may  result  from 
the  stiffness  variation  of  the  clamped  workpiece.  A  work-chuck 
assembly  with  three  clamping  jaws  is  simulated  by  using  a  time 
varying  stiffness  ratio 


(k  /  k  )(1  +  e  sin(3Qt  +  Y  )) 
c  m 


e  -  .2,  Y  -  .8 


Fig.  5.5  shows  that  the  self-tuning  controller  can  compensate  for 
this  kind  of  error  source  as  nicely  as  it  does  for  the  former 
periodic  disturbance.  When  the  phase  angle  Y  is  varied  to  any  other 
value,  the  simulation  results  remain  very  similar.  This  is  an 
advantage  of  using  STC  because,  in  practice,  different  diameters  of 
the  workpiece  have  different  phase  shifts  of  the  jaws  [81], 

The  polar  plot  of  the  uncontrolled  geometrical  error  has  six 
lobes  instead  of  three  in  Fig.  5.5.  This  is  due  to  the  loss  of 
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Figure  5.5  GAC  System  With  Periodic  Disturbance  From 
Stiffness  Variation  of  Chucked  Workpiece 
(Simulation  1-3) 


contact  which  occurs  in  the  last  two  revolutions.  Same  nonlinear 
effect  can  also  be  seen  in  Fig.  5.3,  which  has  four  lobes  rather  than 
two. 

To  see  if  the  self-tuning  controller  can  compensate  for  two  or 
more  periodic  disturbances ,  two  sinusoidal  functions 

.004  sin(flt  +  Y)  +  .002  sin(8.50t) 

are  simulated  as  shown  in  Fig.  5.6.  It  shows  that  the  geometrical 
error  is  reduced  somewhat  but  not  satisfactorily  enough,  if  the  same 
self-tuning  model  and  model  order  are  used.  The  time  varying 
estimates,  together  with  the  coarse  control  accuracy,  indicate  that 
the  assumed  model  order  is  insufficient.  As  we  have  shown  that  the 
cancellation  of  the  periodic  disturbance  is  due  to  the  estimated 
control  model  having  the  same  dynamics  as  that  of  the  disturbance 
function,  therefore,  the  model  order  must  increase  to  four  when  there 
are  two  sinusoidal  disturbances.  Under  the  same  simulation 
conditions  except  that 

n  «  4 

and 

X  -  .99  -  .5  /  t 

Fig.  5.7  shows  that  the  geometrical  error  significantly  reduces  and 
the  parameter  estimation  converges  fast  again.  Besides  these,  the 
simulation  shows  that  the  frequency  of  the  periodic  disturbance  need 
not  be  integer  multiples  of  the  angular  speed  of  the  spindle. 
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It  is  sometimes  convenient  and  desirable  to  offset  the  sensing 
area  from  the  tool  tip  during  the  machining  operation.  The  system 
time  delay  will  be  greater  than  one  in  this  circumstance.  When  k  - 
3,  simulation  shows  that  the  self-tuner  we  used  is  not  fast  enough  to 
minimize  the  transient  response  due  to  the  abrupt  change  of  slope  of 
the  depth-of-cut  in  the  last  but  one  revolution,  and  is  unstable  in 
the  last  revolution  where  loss-of-contact  occurs  most  often.  A 
rational  solution  is  to  alleviate  the  change  of  slope  by  introducing 
a  perturbation  input  as  shown  in  Fig.  5.1c,  where 

6  t-(N-2p) 

A  ( t)  *  **2^'  [!•  +  cos  11 - p - “] ,  if  N-3p  <  t  <  N-p 


0,  otherwise 


(3.3) 


The  perturbation  provides  excitation  for  activating  the  tuning  of  the 
controller  parameters  to  cope  up  with  the  input  changes.  Fig.  5.8 
shows  that  before  applying  perturbation  A(t)  the  controlled  response 
is  still  perfectly  regulated  despite  the  large  time  del^y.  The 
perturbation,  then,  starts  at  the  time  instant  N-3p.  The  transient 
response  is  significantly  reduced  and  satisfactory  out-of-roundness , 
0.00125,  is  thus  obtained.  Since  the  perturbation  tends  to  slow  down 
the  infeed,  the  tool  should  stay  at  the  final  position  rather  than 
retract  at  the  last  revolution. 


The  simulation  of  the  above  deterministic  systems  clearly 
illustrates  how  the  self-tuner  behaves  to  achieve  the  optimal 
performance.  Next  we  shall  find  how  the  same  controller  reacts  in 
the  presence  of  stochastic  disturbances. 
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5.2.2  Stochastic  Systems 


By  adding  the  following  stochastic  disturbance 


W  =  e  +  2.0e  +  1.31e  „  +  .28e 

t  t  t-1  t-2  t-3 

e  :  Gaussian  white  noise  with  zero  mean  and  rms  .002 
t 

to  the  same  machining  system,  Fig.  5.9  shows  that,  compared  to  the 
uncontrolled  result,  the  controlled  dimensional  error  is  minimized 
significantly  and  the  out-of-roundness  and  finish  are  also  reduced  by 
half.  The  controlled  finish,  actually,  is  very  close  to  the  rms  of 
e  ,  the  minimum  that  an  optimal  controller  can  achieve.  The  periodic 
disturbance  is  smeared  by  the  stochastic  disturbance  as  can  be  seen 
in  the  plots  of  the  time  response.  The  estimated  parameters, 
therefore,  do  not  reflect  the  dynamics  of  the  sinusoidal  disturbance 
any  more.  Instead,  they  constitute  a  suboptimal  controller  to 

minimize  the  variance  and  the  bias  of  the  output.  Increase  of  the 
model  order  to  four,  full  order  of  the  simulated  system,  will 
slightly  improve  the  geometrical  accuracy.  Nevertheless  the  true 
parameters  of  the  system  can  still  not  be  identified.  This  is  the 
general  characteristic  of  the  implicit  type  of  self-tuning  control 
when  the  desired  response  is  not  spectrally  rich  enough.  Besides, 
the  insignificant  improvement  in  accuracy  by  the  full  order  means 
that  the  implicit  type  of  controller  is  rather  robust  for  a  mode  of 
lower  order.  For  higher  level  of  disturbance,  Table  5.2  shows  that 
the  self-tuner  maintains  near  optimal  accuracy. 
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Table  5.2  Simulation  Results  of  Different  Disturbance 
Levels  and  System  Delays 


simulation! 

DIMENS 

OUT-O  F- 

FINISH 

condition  | 

ERROR 

ROUND 

[RMS] 

1-7  1 

0.00020 

0.01074 

0.00242 

1 

(-.00162) 

(0.02097) 

(0.00552) 

1-8* 

0.00070 

0.02217 

0.00469 

|  (-.00048) 

(0.04205) 

(0.00955) 

1-9* 

0.00226 

0.04411 

0.00879 

(0.00360) 

(0.07682) 

(0.01801) 

1-10** 

0.00161 

0.02009 

0.00444 

(-.00162) 

(0.02097) 

(0.00552) 

** 

1-11 

0.00122 

0.01877 

0.00469 

(-.00162) 

(0.02097) 

(0.00552) 

(data)  -  result  of  uncontrolled  system 

*  -  high  level  of  disturbance  (see  Table  5.1) 

**  -  large  time  delay  (see  Table  5.1) 
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When  the  system  delay  k  is  greater  than  one,  as  we  have  learned, 
perturbation  in  the  last  few  revolutions  of  machining  is  necessary  to 
alleviate  the  input  changes  and  make  the  controller  stable.  This  is 
particularly  true  in  a  stochastic  system.  Only  with  the 
perturbation,  (5.3),  we  can  obtain  the  geometrical  accuracies  as 
shown  in  Table  5.2  for  k«2  and  3.  The  results,  though  quite  poor 
compared  with  those  for  k*l,  are  still  near  optimal,  because  the 
theoretical  rms  for  k*2  is  0.00448,  which  coincides  with  the 
simulated  results.  The  change  of  the  output  variance  due  to  the 
increase  in  time  delay  may  not  be  so  large  as  was  obtained  through 
simulation.  It  depends  on  the  actual  dynamics  of  the  stochastic 
disturbance  and  the  system  in  real  operation.  A  good  practice  is  to 
reduce  the  system  delay  for  measurement  as  much  as  possible. 

The  important  phenomenon  of  the  machining  process,  chatter,  can 
be  simulated  by  increasing  the  stiffness  ratio.  It  is  found  by 
simulation  that  the  uncontrolled  system  begins  to  be  unstable  as  the 
stiffness  ratio  increases  to  .15,  which  is  close  to  the  theoretical 
value  0.105.  Fig.  5.10  shows  that,  when  the  ratio  is  0.20,  the 
uncontrolled  system  becomes  more  unstable  and  produces  regenerative 
chatter.  And,  as  the  theory  in  Chapter  3  predicted,  the  chatter 
frequency  is  10.8  Hz  which  is  slightly  larger  than  the  natural 
frequency,  9.55  Hz.  The  controlled  system,  however,  remains  stable 
with  accuracy  as  satisfactory  as  that  shown  in  Fig.  5.9.  The 
stiffness  ratio  can  even  be  raised  to  6.0,  as  shown  in  Fig.  5.11,  and 
the  controlled  system  maintains  good  geometrical  accuracy. 
Therefore,  the  increase  in  stability  is  at  least  by  a  factor  of  40 


M 


(*6.0/0.15)  in  Che  simulation.  It  is  comparable  with  the  results  of 
Comstock  et  al.  [29]  and  of  Nactigal  and  Maddux  [76].  According  to 
the  analysis  in  Chapter  4,  since  the  simulated  system  has  a  stable 
polynomial  B(z  *)  in  (4.1),  the  controlled  system  is  stable  no  matter 
how  large  the  stiffness  ratio  is.  This  is  similar  to  what  is  claimed 
by  Mitchell  [66]  that  a  factor  of  1000  increase  in  stability  is 
attainable.  Experimentally,  however,  the  stiffness  ratio  is  found  to 
be  of  the  order  of  0.5  to  .02  and  damping  ratio  to  be  about  0.25  to 
0.01  [104].  Therefore,  it  is  more  appropriate  to  say  that  an 
accurate  and  chatter-free  machining  can  be  obtained  by  means  of  the 
self-tuning  controller. 

To  see  the  necessity  of  incorporating  the  physical  insight  into 

the  self-tuning  model,  the  general  model  in  Chapter  3  is  tested  and 

compared.  A  constant  but  unknown  term  is  added  in  the  model  to 

account  for  the  static  deflection  resulting  from  the  cutting  load. 

With  the  same  order  of  model  and  other  simulation  conditions  as  those 

of  Fig.  5.10,  the  general  self-tuner  also  results  in  a  stable  system 

as  shown  in  Fig.  5.12.  The  controlled  accuracies,  however,  are  much 

worse  than  those  of  Fig.  5.11,  which  uses  the  proposed  model  and 

operates  in  a  more  adverse  condition,  k  /k  *6.0.  It  is  found  that,  as 

cm 

the  stiffness  goes  higher,  the  controlled  inaccuracies  of  the  general 
self-tuner  become  unreasonably  large  and  the  system  can  be  considered 
as  being  unstable.  Hence,  careful  modeling  of  the  process  to  be 
controlled  is  very  important  to  enhance  the  effectiveness  of  self¬ 
tuning  control. 
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5.3  Simulation  With  Multiprobe  Measurement 


The  second  part  of  the  simulation  is  to  test  the  GAC  system  with 
multiprobe  measurement.  The  interesting  properties  of  the  multiprobe 
measurement,  as  we  have  analyzed  in  Chapter  4,  can  be  verified  by 
simulation  with  deterministic  disturbance  only.  Selecting  proper 
value  of  the  detuning  factor  is  very  important  to  achieve  geometrical 
accuracies  as  good  as  those  of  the  former  simulation.  With  the 
stochastic  disturbance  added  in,  the  performance  of  noise  rejection 
and  Increase  in  stability  can  be  evaluated. 


5.3.1  Configuration  of  Probes  for  Simulation 


The  configuration  of  the  multiprobe  measurement,  as  shown  in 
Fig.  4.3,  is  as  follows: 


o  o 

a  -  90  ,  S  -  30 


and  thus,  (4.21)  gives 


,  -1  2  -5  1  -20 

P,  (z  )  -  -  -  Z  +  -  z 

d  &  J3 


(5.4) 


which  has  the  following  characteristic  roots: 


0.933  ej9,  6  -  ±25.55°,  ±46.45°,  ±97.55°, 


±118.45°,  or  ±169.55° 


10  0  0  0  o  o 

1.000  e  ,  0  «  ±6  ,  ±66  ,  ±78  ,  ±138  ,  or  ±130 


This  arrangement  is  favorable  for  the  self-tuning  control  because  all 

the  zeros  of  P,  (z  S  lie  inside  or  on  the  unit  circle.  However,  from 
d 

*  * 

(4.30),  it  suppresses  harmonics  with  frequencies  m  0  ,  where  m 

equals  to 


m*  ’  m°d[g^d(m^,mg)  1  l*  P1’  “a  “  15,  me  -  5,  p  -  60 


mod[ 12k  ±  1,  60],  k  -  0,  1,  2,... 


1,  11,  13,... 


We  shall  find  what  the  controlled  geometrical  errors  would  be,  if  the 
periodic  disturbance  happens  to  be  suppressed  by  the  multiprobe 
measurement.  Solutions  will  be  made  to  overcome  the  control 
deficiency  in  this  type  of  situation. 


5.3.2  Test  of  Detuning  Factor 


As  stated  in  Chapter  4,  p, (z  )  has  at  least  one  zero  on  the 

d 

unit  circle;  therefore,  a  detuning  factor,  0  <  p  <  1 ,  is  needed  in 
(4.35)  for  computing  the  control  input.  Without  using  the  detuning 
factor,  i.e.,  p-1,  Fig.  5.13  shows  unsatisfactory  control  results 

under  the  same  conditions  of  simulation  as  shown  in  Fig.  5.3.  We  see 
that  the  uncontrolled  response  is  replaced  by  the  controlled  response 
which  has  compound  sinusoidal  wave  behavior  during  cutting.  This 
compound  wave,  through  spectral  analysis,  consists  of  harmonics  with 
frequencies  ft,  lift,  and  13ft,  which  are  the  components  invisible  to 
the  multiprobe  measurement.  When  the  detuning  factor  is  added  in. 


P*.95,  Fig.  5.14  shows  that  this  compound  wave  dies  out 
asymptotically,  and  satisfactory  geometrical  accuracy  is  obtained 
thereafter.  The  controlled  response  also  shows  that,  due  to  the 
detuning  factor,  the  controller  acts  less  optimally.  The  specified 
infeed  is  not  tracked  as  perfectly  as  that  in  Fig.  5.3  during  the 
operation.  Therefore,  instead  of  retracting  the  tool,  the  tool  is 
commanded  to  stay  at  the  specified  position  in  the  last  revolution  of 
cutting.  This  gives  the  self-tuner  more  opportunities  to  reduce 
excess  geometrical  errors. 

5.3.3  Test  of  Invisible  Harmonics 

For  the  case  of  suppression  of  harmonics  in  multiprobe 
measurement,  Fig.  5.15  shows  that  the  geometrical  errors  resulting 
from  periodic  disturbance  with  frequency  ft  can  not  be  removed  from 
the  GAC  system.  This  result  is  quite  understandable  because  the 
multiprobe  measurement  is  unable  to  detect  any  errors  having  this 
frequency.  The  simulation,  however,  shows  that  the  higher  harmonics 
in  the  uncontrolled  polar  plot,  which  result  from  the  nonlinearity  of 
the  system,  can  be  compensated  for  by  the  controller.  The  geometry 
of  the  controlled  machined  part  is  a  near  perfect  but  eccentric 
circle.  The  indicated  out-of-roundness  actually  represents  twice  the 
eccentricity  of  the  machined  part.  For  products  whose  roundness  is 
much  more  critical  than  the  eccentricity,  the  GAC  system  with 
multiprobe  measurement  is  still  very  effective.  Otherwise,  the  best 
way  is  to  balance  this  periodic  disturbance  before  the  cutting 


operation  takes  place. 


aw 


( UNCONTROLLED ) 
DIHENS  ERROR  -.00123 
OUT-OF-ROUND  .00861 


(CONTROLLED) 
OWENS  ERROR  ■ 
OUT-OF-ROUNO 
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0  90  180  270  360  *50  540  630  720  810  900 

Figure  5.14  The  Effect  of  Detuning  Factor  in  GAC  System 
With  Multiprobe  Measurement 
(Simulation  1 1—  2 ) 
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The  suppression  of  the  fundamental  harmonic  is  common  to  all 
configurations  of  the  multiprobe  measurement.  The  next  higher 
suppressed  harmonic  of  this  particular  configuration  is  of  frequency 
1 lft .  Conceivably,  Fig.  5.16  is  obtained  by  showing  that  the 
reduction  of  the  periodic  disturbance  is  very  small.  To  overcome 
this  deficiency,  the  configuration  is  changed  to 

o  o 
a  -  132  ,  6  -  30 


that  is 


-1  -5  -27 

P_,(z  )  -  2.40487  -  z  +  1.61803z 

d 


which  has  unstable  zeros 


1.00288 


e±jl40. 43118° 


1.00444 


±j73.18414° 
e  J 


And,  from  (4.30),  no  harmonics  are  suppressed  by  this  configuration 
except  the  fundamental  one  with  frequency  3. 

-1 

Since  some  of  the  zeros  of  the  new  F^Cz  )  lie  outside  the  unit 

circle,  lower  detuning  factor,  p*.85,  is  required  to  stabilize  the 

control  system  and  result  in  satisfactory  accuracies  as  shown  in  Fig. 

5.17.  However,  the  low  detuning  factor  also  detracts  the  performance 

somewhat.  Other  configurations  can  give  similar  results,  too.  It  is 

-1 

found  by  simulation  that  the  farther  the  unstable  zeros  of  Pd(z  ) 
are  away  from  the  unit  circle,  the  lower  the  detuning  factors  needed 
is,  and  the  less  satisfactory  the  accuracy  can  be.  The  configuration 
=  6  has  all  zeros  of  P,(z  S  on  the  unit  circle,  but  it  is 


with  a 
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Figure  5.16  GAC  System  for  Higher  Harmonic  Invisible  to 
Multiprobe  Measurement  ( Siraulat-ionII-4) 
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Figure  5.17  GAC  System  With  Modified  Configuration  of 
Multiprobes  (Simulation  II-5) 


very  possible  to  miss  the  detection  and  the  control  of  the 

disturbances  having  frequencies  not  equal  to  integral  multiple  of 

frequency  ft.  The  determination  of  the  probe  configuration  depends 

on  the  convenience  of  the  spatial  arrangement  of  the  probes,  the 

prior  knowledge  of  the  disturbance  frequencies,  and  the  locations  of 

the  zeros  of  P  (z  *)• 
d 

5.3.4  Test  of  Stochastic  System  and  Stability 

Here,  we  return  to  the  original  probe  configuration 

o  o 

corresponding  to  O“90  and  &»30  .  Under  the  same  conditions  of 
simulation  as  those  used  in  Fig.  5.9,  good  control  accuracies  are 
obtained  as  shown  in  Fig.  5.18  in  the  presence  of  stochastic 
disturbance.  The  selected  detuning  factor,  p«0.88,  is  smaller  than 
that  used  in  Fig.  5.14,  equal  to  0.95.  In  general,  a  stochastic 
control  system  requires  a  lower  detuning  factor  t  an  a  deterministic 
system  does.  The  forgetting  factor,  however,  changes  in  the  opposite 

direction  because  the  self-tuning  controller  needs  more  data  for 

variance  reduction  in  the  stochastic  cases.  Also,  it  is  found  by 
simulation  that  a  large  forgetting  factor  is  good  for  a  low  detuning 
factor.  Therefore,  the  forgetting  factor 

X  -  1.0  -  .5  /  t 

is  used  in  Figs.  5.17  and  5.18  and  the  following  simulation. 

As  the  stiffness  ratio  increases  to  0.2,  the  controlled  system 
remains  stable  and  accurate,  again.  Comparing  the  controlled 

responses  of  Fig.  5.19  to  those  of  Fig.  5.18,  we  see  that,  in  the 
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Figure  5.19  GAC  System  With  Multiprobe  Measurement  for 
Unstable  Machining  Process 
(Simulations  11-7) 
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middle  of  the  operation,  the  static  error  in  Fig.  5.19  is 
approximately  twice  that  of  Fig.  5.18  with  the  same  detuning  factor. 
And,  as  shown  in  Fig.  5.20,  this  static  error  increases 
proportionally  when  the  stiffness  ratio  goes  up  to  2.0.  The  system  is 
still  stably  controlled.  However,  in  the  stage  of  dwell  cutting,  the 
static  error  decreases  more  slowly  than  the  specified  depth-of-cut . 
Consequently,  increase  of  the  duration  of  dwell  cutting  will  result 
in  better  geometrical  accuracy.  In  Fig.  5.20,  an  extra  revolution  of 
dwell  cutting  is  added  and  satisfactory  accuracies  are  obtained  as 
shown.  The  large  dimensional  excess  suggests  that  another  revolution 
of  dwell  cutting  can  be  beneficial. 

The  increase  in  stability,  found  by  simulation,  is  by  about  a 
factor  of  20.  The  detuning  factor  is  the  main  reason  for  the  limited 
closed  loop  stability.  However,  the  improvement  in  performance  is 
quite  significant. 

5.4  Summary  of  Discussions 

In  the  first  part  of  the  simulation,  where  the  measurement 
problem  is  not  considered,  we  observed  the  following: 

(I. a)  Deterministic  disturbances  can  be  fully  compensated  for  by  the 
self-tuning  controller  irrespective  of  the  amplitudes,  the  frequency 
contents,  and  the  types  of  the  disturbances. 

(I.b)  A  time  varying  forgetting  factor  is  very  usetul  not  only  to 
speed  up  the  convergence  at  the  startup  stage,  but  also  to  enhance 
the  adaptability  of  the  controller  for  the  environmental  changes  or 
modeling  errors. 


(I.c)  Self-tuning  controller  gives  perfect  tracking  during  the 
operation.  However,  if  the  system  time  delay  is  greater  than  one, 
the  controller  may  not  be  able  to  follow  the  input  change  as  quickly 
toward  the  end  of  the  cutting  operation.  Perturbation  can  be  added 
before  and  during  the  input  change.  Transient  response  will  thus  be 
reduced. 

(I.d)  In  the  presence  of  various  levels  of  deterministic  disturbance 
and  stochastic  disturbance,  the  implicit  type  of  self-tuning  control 
is  successful  in  accomplishing  near  optimal  performance  with  lower 
order  model.  The  controller  we  used  is  robust. 

(I.e)  The  increase  in  stability  can  be  by  a  factor  of  at  least  40  for 
the  self-tuning  controller.  In  other  words,  the  self-tuned  system  is 
chatter  free. 

(I.f)  A  self-tuning  model  derived  with  physical  insight  into  the 
machining  process  can  give  much  better  performance  than  the  general 
self-tuner. 

In  the  second  part,  multiprobe  is  used  to  tackle  the  measurement 
problem  of  the  GAC  system.  By  simulation,  we  have  the  following 
important  findings: 

(II. a)  A  detuning  factor  is  crucial  for  the  success  of  the  GAC  system 
with  multiprobe  measurement.  Normally,  it  has  a  positive  value 
smaller  than  one. 

(II. b)  Hie  multiprobe  measurement  always  suppresses  certain 
harmonics.  Those  disturbances,  thus,  cannot  be  reduced.  For  the 
disturbance  with  the  same  frequency  as  the  spindle  rotation,  the  GAC 
system  can  not  do  much  except  removing  disturbance  induced  by  the 


nonlinearity  of  the  system  and  generating  a  perfect  round  but 
eccentric  part.  For  higher  harmonics,  the  uncontrollability  can  be 
avoided  by  selecting  the  probe  configuration  properly. 

(II. c)  The  probe  configuration  determines  the  zeros  of  Pj(z  )•  The 
farther  the  unstable  zeros  are  away  from  the  unit  circle,  the  smaller 
the  detuning  factor  used  should  be,  and  the  less  satisfactory  the 


geometrical  accuracies  obtained  are. 

(II. d)  Stochastic  systems  need  a  lower  detuning  factor  than 


j  „  _ ~  _ _  _ j _ r _ _ _ _ i _ _ ^  - - 4.1 _ 


CHAPTER  6 


A  COMPARISON  OF  STC  AND  FCC 

As  stated  in  Chapters  1  and  4,  the  stochastic  properties  of 
machining  processes  should  be  considered  while  modeling  a  GAC  system 
for  removing  the  nonrepeatable  machining  errors  [47,48,77,7b], 
However,  research  engineers  have  just  started  to  accomplish  this 
actively  in  recent  years.  Rao  and  Wu  [84]  proposed  "compensatory 
control"  for  roundness  control  in  cylindrical  grinding.  Of  late  this 
technique  was  expanded  and  called  "forecasting  compensatory  control" 
(FCC)  [54,55,68,69]. 

Basically,  FCC  is  a  GAC  method  which  uses  in-process 
measurement.  The  machining  process  is  modeled  by  a  stochastic 
process,  AR(m) ,  in  FCC.  The  compensation  is  obtained  by  giving  an 
input  which  cancels  the  error  predicted  by  the  stochastic  model.  The 
modeling  strategy  of  FCC  ignores  the  relationship  between  the  inputs 
and  the  outputs  of  the  machining  system.  Also,  it  lacks  a  systematic 
approach  for  designing  the  compensatory  controller.  Thus,  the 
performance  of  the  compensated  system  can  hardly  be  assured. 

Realizing  this,  we  propose  the  STC  method  in  Chapter  4  for 
cylindrical  machining  processes.  This  technique  adopts  a  more 
rigorous  modeling  and  control  strategy.  The  stochastic  model  is 
deduced  from  the  machine  tool  dynamics  which  has  been  discussed  in 


2 


Chapter  2.  Therefore,  the  self-tuning  controller  can  be  designed  and 
analyzed  systematically  for  a  better  control  performance. 

The  objective  of  this  chapter  is  to  compare  the  performance  of 
these  two  control  strategies  by  theoretical  investigation,  digital 
simulation,  and  analog  test.  In  theoretical  investigation,  the 
theoretical  control  performances  of  both  control  techniques  are 
examined  in  the  aspects  of  command  tracting  and  disturbance 
rejection.  In  digital  simulation,  a  metal  cutting  process  is 
simulated  to  test  the  machining  accuracy  and  stability  obtainable 
from  both  strategies.  The  analog  test  is  to  compare  the  control 
performances  of  both  strategies  which  are  subject  to  the  constraints 
of  finite  word  length  and  sampling  in  microprocessor  based  control 
systems . 


6.  i  Theoretical  Investigation 

For  simplicity,  a  general  single  input  single  output  (SISO) 

system,  as  shown  in  Fig.  6.1,  is  used  for  this  comparison.  All  the 

variables  and  dynamics  are  presented  in  discrete  time  domain.  The 

dynamics  of  the  process,  G  ,  is  assumed  to  have  the  following 

P 

transfer  function 


where 


G  =  B(z_1 )  /  A(z_1) 
P 


. .  -1 .  ,  -1  -n 

A (z  )  =  1  a  z  +  ...  +  a  z 
1  n 


(6.1) 


-lx  -k  „  -1 

B(z  )  =  z  B(z  ) 


A.--  ■■■  v 


- ;  v- 


and  z  is  a  backward  shift  operator.  The  coefficients,  a^'s  and 

b^'s,  are  related  to  the  dynamic  characteristics  of  the  process,  and 

they  are  considered  to  be  unknown  in  adaptive  control.  The  order  n 

and  the  time  delay  k  are,  however,  assumed  to  be  known  here.  The 

deterministic  disturbance  z  and  the  stochastic  disturbance  w  are 

t  t 

unmeasurable,  but  bounded.  The  observed  output  y  is  fed  to  the 

t 

controller  G  ,  where  the  control  input  u  is  generated  by  the 
c  t 

specified  control  strategy  or  algorithm. 

STC 

The  design  of  the  STC  controller  is  based  on  the  certainty 
equivalent  principle.  Therefore,  we  shall  first  discuss  the  design 
of  the  optimal  controller  as  a  convenient  way  for  investigating  the 
theoretical  control  performance  of  the  STC  sysetem.  A  theorem  for 
the  STC  controller,  then,  will  be  presented  to  substantiate  the 
f indings. 

The  output  of  the  system  as  shown  in  Fig.  6.1  is  governed  by 

y  =Gu  +z  +  w  (6.2) 

t  p  t  t  t 

To  simplify  the  comparison,  the  stochastic  disturbance  w  can  be 

t 

omitted  presently  without  loss  of  generality.  Hence, 

y=Gu+z  (6.3) 

t  p  t  t 


Since  most  of  the  stationary  deterministic  disturbances  are  constant 


or  periodic  in  the  time  domain,  one  can  always  find  an  operator  of 
the  form 


such  that 


-1  -1  -s 

P(z  )  =  p  +  p  z  +  .. .  +  p  z 
0  1  s 


P( z  *  )z  =0 


(6.4) 


By  substituting  (6.1)  into  (6.3),  it  gives 


B'(z_1) 

y  *»  - : -  u  ,  +  z 

c  . .  -1.  t-k  t 


(6.5) 


A(  z  ) 


Suppose  that  the  output  y  can  be  predicted  bv  the  linear 

t+k 


combination  of  the  observed  outputs  and  inputs,  that  is, 


-1  -1 

y  ,  =  F(  z  )y  +  G(  z  )u 

-’t+k  }  t  t 


(6.6) 


where 


F(z  S  =  fQ  +  f2z  1  +  ... 


-1  -1 

G(z  )  =  gQ  +  gjZ  +  ... 


Thus,  for  the  desired  output  y  the  control  input  should  be 


ut  “  ,  -L  yt+k 


F(z 


G(  z  ) 


G(z  *) 


-1  s  yt 


(6.7) 


Substituting  (6.7)  into  (6.5)  gives 


t+k  AG  +  BF  t+k.  AG  +  BF  t+k 


where  the  operator  z  is  omitted  for  simplicity  in  the  expression. 

* 

In  order  to  achieve  y  ,,  *  y_,,  ,  F  and  G  must  satisfy 

t+k  t+k 

B'  -  AG  +  BF  =  AG  +  z  B'F  (6.9) 

Gz  -  0  (6.10) 

t+k 

From  (6.4),  if  P  is  a  factor  of  G,  the  condition  (6.10),  cancellation 
of  the  deterministic  disturbance,  will  be  satisfied.  Also,  the 

condition  (6.9)  requires  that  there  be  a  factor  B'(z  *)  in  G(z  *)• 
Therefore, 


G  -  EB'P 


(6.11) 


where 


E  =  e  +  e  z  *  +  . . . 

0  1 

Substituting  (6.11)  into  (6.9)  gives 

1  -  APE  +  z~kF  (6.12) 

There  exists  a  unique  solution  for  E  and  F  in  (6.12),  if  the  order  of 
E  equals  to  k-1,  and  the  order  of  F  equals  to  n+s-1. 

The  above  derivations  demonstrate  that  the  system  can  be 

* 

optimally  controlled,  y^  *  y^,  if  the  controller  (6.7),  with  F  and  G 
satisfying  (6.11)  and  (6.12),  is  used.  In  STC,  where  the  polynomials 


F  and  G  have  to  be  estimated  through  a  certain  adaption  mechanism, 
the  same  optimality  can  also  be  achieved.  This  can  be  seen  from  the 
following  theorem  [39]: 


Theorem:  If  the  system  (6.5)  satisfies 

(i)  all  the  zeros  of  BP  lie  inside  or  on  the  unit  circle. 

(ii)  all  the  zeros  of  transfer  function  B/A  lie  strictly 
inside  the  unit  circle. 

(iii)  the  zeros  of  the  polynomial  BP  on  the  unit  circle  have 
a  Jordon  block  of  size  one. 

the  self-tuning  controller 


yt+k  =  Fyt  +Gut  (6*13) 


e-s  +  S  $  (y  —  <}>  0  )  (6.14) 

t  t-1  t-1  t-k  t  t-k  t-1 


S  *  <b  s 

,  _  .  t-2VkVk  t-2  . 

st-i  e-2  ■  ~  (6*15) 

1  +  ^t-kSt-2^t-k 

where 


9t  ■  tro,r  1  •  —  ’®o’ej  ■ — 1 


(t  ■  [VVi . Wir 


will  lead  to  a  stable  closed  loop  system,  and  the  output  will 
* 

converge  to  yt  asymptotically. 

The  inclusion  of  the  polynomial  P  in  the  above  assumptions, 
which  results  from  the  existence  of  the  deterministic  disturbance  z  . 


can  be  taken  away  if  P  has  no  multiple  roots  on  the  unit  circle.  For 
the  system  having  the  stochastic  disturbance  in  addition  to  the 


deterministic  disturbance,  the  above  self-tuning  controller  can  be 

extended  by  incorporating  extra  parameters  in  6  to  account  for  the 

t 

dynamics  of  the  stochastic  disturbance.  Again,  the  disturbances  can 
be  reduced  and  the  specified  output  will  be  tracked  [39]. 

An  important  point  of  the  theorem  is  that  the  estimated 
polynomials  F  and  G  are  not  necessary  to  satisfy  the  relationship  as 
(6.12)  for  the  optimal  performance  of  the  output.  This  usually 
happens  in  the  applications  of  STC,  particularly  for  the  above  type 
of  self-tuning  controller  which  estimates  directly  the  parameters  of 
the  controller  rather  than  those  of  the  process  to  be  controlled. 
The  STC  system  will  find  its  own  parameters  to  satisfy  the  control 
objective.  This  can  be  seen  from  the  simulation  results  in  the  next 
section. 

FCC 


In  FCC,  it  is  assumed  that  the  output  of  an  uncontrolled 
machining  system,  y^ ,  can  be  fitted  by  an  autoregressive  model,  AR(ra) 
[83] 


$'y  “  e 

t  t 


-1  -m 

"  1  -  '  z  -  ...-$' z 

1  ra 


This  model  can  be  derived  from  (6.2)  by  letting  u  “0  and  converting 

t 

the  disturbances  into  an  autoregressive  process  with  white  noise  e  . 


By  the  following  factorization 


1  “  t'E  +  z 


-k. 


-1  -k+1 

E  ■  e,.  +  e,  z  +  . . .  +  e.  .  z 
0  1  k-1 


-1 

4  +  A  Z  + 


■m+ 1 


m-1 


the  output  y  can  be  predicted  by  the  linear  combination  of  its 
previous  outputs  at  time  t,  t-1,...,  t-m,  that  is, 


't+k 


^t+k  + 


Ee 


t+k 


>t+k  *  %  (6'16) 

* 

For  the  desired  output  y  ,  ,  the  control  input  of  FCC  is  [84,54,55, 

t+k 

68,69] 


u 

t 


yt+k  ~  y t+k  *  yt+k 


4>y 


(6.17) 


For  the  system  with  deterministic  disturbance  only,  substituting 
(6.17)  into  (6.5)  gives 


t+k 


A  +  B'  y t+k  +  A  +  B' 


t+k 


(6.18) 


Equation  (6.18)  shows  that  it  is  not  possible  to  realize  y  -  y 

t+k  t+k 

or  to  cancel  z  ,  provided  that  the  control  law  (6.17)  is  used.  The 
t+k 

reason  for  the  inability  to  fully  track  the  specified  output  or 
eliminate  the  disturbance  is  that  FCC  does  not  consider  the  dynamic 


29 


relationship  between  the  system  output  and  the  compensation  input. 
This  is  not  the  case  for  the  in-process  control  of  the  machining 
system. 


* 

When  y  »*0, 
t+k 


from  (6.18), 


it  gives 


A 

yt+k  A  +  B'  2 t+k 


(6.19) 


The  above  closed  loop  system  is  the  same  as  a  conventional  control 
system  with  a  unit  negative  feedback.  The  simulation  in  the  next 
section  will  show  this  characteristic. 

6.2  Digital  Simulation 


In  digital  simulation,  the  plunge  cutting  process  described  in 
Chapters  4  and  5  is  simulated  and  both  techniques  are  tested  on  a  CDC 
6000  computer  system. 

The  model  used  in  STC,  from  (4.1),  is 


where 


A(y  -  y  +  <5  )  **  B(u  -  y  +  6  )  +  Ce  (6.20) 

t  t-p  t  t-k  t-p  t  t 


-1  -n 

A  *  1  +  a  z  +  . . .  +  a  z 
1  n 


B  *  b.  +  b.z  ^  +  ...  +  b  z  n 
0  1  n 


-1  -n 

C  *  1  +  c  z  +...+cz 
1  n 


The  STC  algorithm,  from  (4.11)  and  (4.12),  is  to  update  the 


parameters  of  F,G,  and  C  in  the  following  model 


y  *  F(y  -  y  +6  )  +  G(u  -  y  +  6  )  + 

t  t-k  t-p-k  t-k  t-k  t-p  t 


C(y  -  6  )  +  e 

t-p  t  t 


(6.21) 


and  to  generate  the  control  input  by  the  following  control  law 


u 

t 


(-1/g  )[F(y  -  y  +  6  )  +  (G  -  g  )(u 

0  t  t-p  t  0  t 


y  ,  +  6  ) 

t-p+k  t+k 
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In  designing  the  FCC  controller,  the  model 
estimation  is 


(6.22) 


for  parameter 


y  “  $  y  +  $  y  +  ...  +  ®  y  +  $  +  c  (6.23) 

yt  Vt-k  Vt-k-1  Vrt-k-m+1  ro  t 


where  the  estimate  $  is  added  to  account  for  the  nonzero  value  which 

m 


may  exist  in  y  .  And  the  control  law  is 
t 


u 
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Recursive  least  squares  method  is  used  for  updating  the 
parameters  of  both  control  techniques.  Because  of  the  existence  of 
the  deterministic  disturbance  the  order  of  the  model  in  (6.20)  is  not 
less  than  four.  However,  to  test  the  robustness  of  the  STC  system 
for  the  reduced  order  of  the  model,  n*2  is  assumed.  Hence,  there  are 
seven  parameters  to  be  estimated  in  (6.21) 

fg,fj  »§q,8|  ,^2’^1  ’  an<*  ^2*  F°r  Che  Same  comPlexity  in  parameter 
estimation,  therefore,  the  order  ra  in  (6.23)  is  assumed  to  be  six. 


The  simulation  results  depicted  in  Fig.  6.2,  where  the 
stochastic  disturbance  is  set  to  zero,  demonstrate  that  the 


originally  unstable  cutting  system,  k  /k  -0.2,  can  be  stabilized  by 

c  a 

the  self-tuning  controller,  and  the  machining  error  can  be  nearly 
eliminated.  The  polar  plot  represents  the  final  geometrical  error  of 
the  machined  part.  The  control  input  shows  that  the  self-tuning 
controller  automatically  compensates  for  the  unknown  deterministic 
disturbance  and  the  structural  deflection  resulting  from  the  nominal 
cutting  force.  The  estimates  g^  ,  g^  and  g^ ,  which  are  related  to 
the  control  input,  satisfy  (6.4).  It  shows  how  the  self-tuning 
controller  adjusts  its  parameters  to  counter  the  deterministic 
disturbance. 

For  the  stochastic  machining  system  as  shown  in  Fig.  6.3,  STC, 
again,  stabilizes  the  system  and  produces  good  accuracy.  The  error 
left  on  the  machined  part  has  the  root-mean-squares  (rms)  value 
approximating  to  that  of  e^  in  Chapter  5,  0.002,  which  is  the  minimum 
value  we  can  get. 

When  the  structural  compliance  is  doubled,  k  /k  »0.4,  the 

c  m 

uncontrolled  cutting  system  becomes  more  unstable  as  shown  in  Fig. 
6.4.  The  STC  system  still  shows  good  stability  and  accuracy. 

As  mentioned  in  the  theoretical  investigation,  the  STC  algorithm 
used  in  the  above  simulations  belongs  to  the  direct  or  implicit  type 
of  adaptive  control  [39],  which  tends  to  update  the  estimates  for 
minimizing  the  output  error  rather  than  for  matching  the  true 
parameters.  Hence,  we  can  get  the  optimal  or  nearly  optimal  output 
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despite  the  nonconvergence  of  the  estimates.  Furthermore,  the  linear 
and  reduced-order  controller  provides  a  certain  degree  of  robustness 
by  showing  satisfactory  stability  and  accuracy  in  the  simulated 
system  which  has  a  more  complicated  structure. 

Figs.  6.5  and  6.6  show  the  simulation  results  from  FCC  with  the 
simulation  conditions  corresponding  to  those  of  Figs.  6.2  and  6.3 
respectively.  Although  the  systems  are  stabilized  by  the  forecasting 
compensator,  the  accuracies  of  the  parts  are  much  worse  than  those 
generated  by  the  STC  system.  The  out-of-roundness  shown  in  Fig.  6.5 
is  about  four  times  that  shown  in  Fig.  6.2.  The  control  response  and 
the  control  input  show  that  the  structural  deflection  during  the 
machining  operation  is  only  partly  compensated  for.  It  confirms  the 
conclusion  from  the  theoretical  investigation  that  the  forecasting 
compensator  is  essentially  a  unit  feedback  controller. 

By  considering  the  confidence  interval  of  the  parameter 
estimation,  the  estimates  of  Fig.  6.6  reveal  that  the  FCC  machining 
system,  when  k  =  1,  can  be  modeled  as  a  random  walk  process 

y  *  y  +  e  (6.25) 
t  t-1  t 

Consequently,  the  control  input  is 


u  =  -  y  (6.26) 

t  t 


which  is  the  simple  compensatory  control  (SCC)  by  one  step  delay 
[54],  The  control  inputs  and  the  system  outputs  in  Figs.  6.5  and  b.b 


satisfy  approximately  the  relationship  given  by  (6.26) 


The  same 


[I, 


CONTROL  INPUT 
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results  can  be  seen  in  the  cutting  test  conducted  by  Moon  et  al. 
[69]. 


For  the  less  rigid  system  as  that  simulated  in  Fig.  6.4,  FCC  is 
unable  to  stabilize  the  system  any  more,  as  shown  in  Fig.  6.7. 
Further  simulation  show  that  the  increase  in  stability  gained  by  STC 
can  be  five  times  that  by  FCC. 

6.3  Analog  Test 

In  the  above  simulation,  we  have  compared  the  control 
performance  of  the  two  control  techniques  by  using  the  machining 
model  of  chatter  control  and  the  experimental  machining  properties. 
Therefore,  similar  control  performance  may  be  expected  in  the  control 
of  the  real  machining  operation.  However,  some  considerations  in 
real  implementation  still  need  to  be  tested  for  assuring  the  actual 
performance.  The  analog  test  is  to  examine  the  control  performance 
of  the  two  control  strategies  subject  to  finite  word-length  and 
sampling  rate  of  the  microprocessor  controller.  The  system  setup  is 
shown  in  Fig.  6.8. 

The  electronic  circuit  for  the  test  is  an  RLC  low  pass  filter 
representing  a  second  order  dynamic  system  as  shown  in  Fig.  6.9.  The 
input  disturbance  z  is  a  sinusoidal  signal  generated  by  a  function 
generator.  Without  the  control  loop,  Fig.  6.10  shows  the  output  of 
the  uncontrolled  system.  The  large  sinusoidal  wave  in  the  figure 
results  from  the  deterministic  disturbance  generated  by  the  function 
generator.  Since  the  circuit  is  not  accurately  made,  a  small  wave 
with  higher  frequency  is  superimposed  to  the  large  wave  in  the 


Electronic  Circuit  for  Analog  Test 


10  Output  of  Uncontrolled  System 


figure.  This  unforeseen  nonlinearity  is  also  a  good  source  to  test 
he  robustness  of  the  control  algorithms.  The  objective  of  the 
control  function  is  to  compensate  for  the  effect  of  the  deterministic 
disturbance  and  to  get  the  minimum  output. 

An  M68000-based  data  acquisition  and  control  system  was  devised 
to  test  both  control  techniques  on  the  electronic  circuit.  The 
digital  controller  has  a  12-bit  analog  to  digital  converter  (ADC)  and 
a  12-bit  digital  to  analog  converter  (DAC).  In  numerical 
computation,  each  16-bit  data  word  is  fixed-point  formatted.  The 
higher  byte  represents  the  integer  part  of  the  data  and  the  lower 
byte  stands  for  the  fraction  part  of  the  data.  The  finite  word- 
length  of  the  processor  and  the  converters  limit  the  maximum 
amplitude  of  the  control  action.  The  maximum  sampling  speed  of  the 
digital  controller  mainly  depends  on  the  number  of  parameters  to  be 
estimated.  About  one  millisecond  is  needed  for  the  control  with  five 
parameters. 

Ou ly  the  output  signal  y  is  needed  by  the  microprocessor 
controller.  The  exact  frequency  and  the  amplitude  of  the  disturbance 
are  unknown  to  the  controller.  However,  the  approximate  range  of  the 
frequency  is  necessary  in  order  to  make  sure  that  the  sampling  rate 
of  the  controller  is  sufficiently  high  for  avoiding  the  aliasing 
problem.  In  this  test,  the  frequency  of  the  disturbance  is  about  10 
Hz,  and  the  sampling  rate  is  100  Hz. 

By  using  the  STC  algorithm  as  shown  in  equations  (6.13-15)  with 
n«*2  and  k«l ,  Fig.  6.11  shows  that  the  variance  of  the  output  is 
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nearly  eliminated,  and  the  disturbance  is  countered  by  the  discrete 
control  input  having  the  same  frequency.  The  output  of  the  FCC 
system,  however,  is  only  slightly  improved  as  shown  in  Fig.  6.12. 
The  forecasting  compensator,  obviously,  is  much  less  effective  than 
the  self-tuning  controller.  Again,  from  Fig.  6.12,  the  forecasting 
compensator  behaves  like  a  unit  feedback  controller. 

6.4  Summary 

In  this  chapter,  we  have  theoretically  shown  that  STC  can 
achieve  the  optimal  control  performance  despite  the  unidentified 
deterministic  and  stochastic  disturbances.  And,  FCC,  which  works 
like  a  unit  feedback  controller,  is  less  effective  than  STC  in 
command  tracking  and  disturbance  rejection. 

The  digital  simulation  shows  that  the  STC  technique  can  not  only 
achieve  the  optimal  response,  but  also  stabilize  the  uncontrolled 
machining  system.  The  FCC  technique,  on  the  other  hand,  produces 
less  satisfactory  output  as  found  in  the  theoretical  investigation. 
And  it  has  less  capability  of  stabilizing  the  uncontrolled  system. 
The  simulation  also  shows  that  the  FCC  system  behaves  as  a  random 
walk  process.  Since  the  experimental  machining  properties  are  used 
in  the  simulation,  similar  results  may  be  expected  in  the  control  of 
the  real  machining  system. 

In  the  analog  test,  the  two  control  techniques  are  tested  and 
compared  by  an  electronic  circuit  with  a  microprocessor  controller. 
The  results  show  that  STC  is  still  more  effective  than  FCC  even  with 


the  constraints  of  the  finite  word-length  and  sampling  rate. 


CHAPTER  7 


CONCLUSIONS  AND  RECOMMENDATIONS 


7.1  Conclusions 

A  new  GAC  system  is  proposed  and  formulated  for  improving  the 
geometrical  accuracies  of  cylindrical  machined  parts.  The 
development  is  accomplished  by  tackling  the  three  major  problems  in 
GAC:  measurement,  modeling,  and  adaptive  control. 

In  measurement,  we  use  Whitehouse's  multiprobe  measurement  to 
solve  the  separation  problem  in  in-process  measurement  of  geometrical 
error.  New  formulations  are  derived  and  analyzed  for  incorporating 
the  multiprobe  measurement  into  the  digital  adaptive  controller. 
This  is  the  first  time  that  multiprobe  measurement  is  used  for 
control  purpose. 

In  modeling,  the  interaction  of  regenerative  cutting  dynamics 
and  structural  dynamics  is  identified  to  be  the  most  prominent 
physics  governing  the  dynamic  behavior  of  cylindrical  machining 
processes.  From  this  we  derive  a  stochastic  model  in  a  discrete 
form  for  the  plunge  cutting  process.  Both  the  stochastic  and  the 
deterministic  disturbances  are  incorporated  in  the  model.  This  model 
is,  though  linear,  very  useful  1  for  designing  a  digital  GAC  system. 


In  adaptive  control,  the  self-tuning  control  theory  is  used  to 
generate  an  adaptive  control  algorithm  for  the  model  we  derived. 
Theoretical  properties  of  the  algorithm  are  analyzed  to  assess  the 
performance  of  the  GAC  system  in  practical  application  and  to  suggest 
the  ways  for  improvement.  Also  derived  is  a  modified  algorithm  for 
the  GAC  system  with  multiprobe  measurement.  A  detuning  factor  is 
introduced  for  stabilizing  the  GAC  system  with  multiprobe 
measurement. 


Simulations  are  carried  out  to  evaluate  the  performance  of  the 
developed  algorithms.  The  results  show  that  the  proposed  GAC  system 
not  only  improves  the  accuracies  considerably  but  also  improves  the 
stability  of  the  machining  system  by  a  factor  of  more  than  forty. 
The  GAC  system  with  multiprobe  has  similar  performance,  except,  as 
the  theoretical  analysis  predicted,  that  certain  harmonic 
disturbances  may  not  be  controllable  for  some  specific  space 
configurations  of  the  probes.  This  problem  can  be  overcome  by 
changing  the  configuration  slightly  and  selecting  an  appropriate 
value  for  the  detuning  factor.  The  increase  in  stability  through  the 
detuning  factor  is  at  the  expense  of  geometrical  accuracies  normally. 
The  simulation  also  shov  that  the  self-tuning  control  algorithms  are 
robust  to  nonlinearities  and  low  model  orders. 


A  comparison  of  STC  and  FCC  shows  that  the  strategies  of 
modeling  and  control  we  used  can  result  in  better  performance  than 
other  existing  strategies.  Also  the  comparison  in  the  analog  test 
shows  that  the  algorithm  of  STC  does  work  well,  as  we  have  predicted, 


in  a  microprocessor  based  control  system. 


7.2  Recommendations 


To  enhance  the  applicability  of  this  new  GAC  system,  the 

following  subjects  are  recommended  for  future  work: 

(1)  A  cutting  test  will  further  substantiate  this  work.  As  we  have 
shown  in  Chapter  6  that  the  microprocessor  controller  is  simple 
to  implement,  the  major  work  will  be  the  setting  up  of  sensors 
and  actuators.  For  internal  grinding  or  turning  processes,  the 
accessibility  of  the  sensors  become  more  critical. 

(2)  In  traverse  cutting  we  need  a  reliable  method  for  determining 
the  overlap  factor  before  the  operation  of  GAC.  In-process 
estimation  of  the  overlap  factor  is  not  practical  because  it 
makes  the  estimation  process  nonlinear. 

(3)  The  strategies  of  modeling  and  control  can  be  modified  and 
applied  to  other  machining  processes,  such  as  milling  [1U2], 
centerless  grinding  [13,85],  EDM,  ECM,  and  laser  machining. 

(A)  Other  sensors,  such  as  force  transducers  and  accelerometers  can 
be  used  in  addition  to  the  displacement  sensor  to  help  the 
measurement  problem  and  to  ensure  the  optimal  performance. 

(5)  The  parameters  estimated  on-line  can  be  used  for  monitoring  the 
conditions  of  the  machining  system  during  machining,  such  as 
tool  wear  or  breakage,  defective  spindle  bearings,  and  other 
malfunctions  of  the  machine  tool. 
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Appendix  1 
Derivation  of  (4.1) 


The  continuous  transfer  function  of  the  tool-workpiece  structure 
can  be  represented  by: 


MT(S)  "  k  2  .  2 

ms  +  2?id  s  +  w 
n  n 


a)  »  natural  frequency  of  the  structure 
n 

C  ■  damping  ratio  of  the  structure 

k  *  stiffness  of  the  structure 
m 

By  taking  the  Z-transformation  with  a  zero-order-hold  [37],  the 
discrete  transfer  function  is  derived  as  follows: 


-1  MT(s ) 

MT  =  (1  -  z  )Z{ — — } 


-1  -2 

j  BjZ  +  &2Z 

r"7T  ir~  “2* 

“  1  ♦  V  °2* 


or  more  generally  [23] , 


1  vvliv" 

MT  =  k  -1  -2 

ra  1  +  <XjZ  +  a^z 


where 


>  AAi 'it  Aj'ji  "jCj  V  T /.  <•.  /.  r.  f.  *.*.■  . 


•  .'»•,.*»  V*  O 


*  -2e  cos  bT,  =  e 


«.gT  3  — 

B,  ■  1  -  e  cos  bT  -  ~  e  sin  bT 
1  b 


-2aT  a  -aT  -aT 

B„  =  e  +  ~  e  sin  bT  -  e  cos  bT 
Z  b 


a  =  5u> 


b  =  “n  V1  '  5 


and  T  is  the  unit  sampling  time  interval. 


The  discrete  form  of  the  regenerative  process  in  Chapter  2  is 


-st 


1  -  pz 


where  the  regenerative  delay  p  is  equal  to  t/T  or  1/(QT).  The 
overlap  factor  is 


0  <5  u  <  1 


The  cutting  process  CP  can  be  represented  by  the  cutting  stiffness, 

i.e.,  CP  *  k  .  The  stochastic  and  the  deterministic  disturbances  of 

c 

Fig.  A. 2  are  given  as 


*  /  .**  m  ■  /*  ***  • 

m.  ■  *  «  "  •  *  ■  **  -»  *  •  *,  •  *  •  m  m 
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w  -  (1  +  Cj'z  +  c2'z  )?c 

Z  =  d  +  d  t  +  Er  sin  (2itiQt  +  6  ) 
t  0  1  i  l 

respectively,  where  5  is  an  uncorrelated  random  process  with  zero 

,  .  2 
mean  and  variance  o^. 

The  GAC  system  is,  thus,  governed  by  the  following  equations 


Y  =  W  +  Z  +  Y  (t)  -  X  +  u  ( A .  1 ) 

t  t  t  s  t  t-k 


-1  -2 
k  0  +  0  z  +  0  z 

X  -  ——2 - 1 - - - (Y  -  pY 

At  k  —i  -2  t  H  t— p 
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m  1  +  V 


+  <*2z 


(A.2) 


By  substituting  (A.2)  into  (A.l)  and  using  the  following  two 
definitions 


y  =  Y  -  Y  (t) 
t  t  s 


6  -  Y  (t)  -  uY  (t-p) 

t  s  s 


we  can  have  the  following  equation 
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where 


(V  +  Yz  +  Yz  )(ut-k  ■  Myt-P  +  V 


+  (b^'  +  b^'z  1  +  b^'z  2)(Wt  +  2^) 


(A. 3) 


6.  =  6  k  /k  ,  i  -  0,1,2 
l  1cm 


a±'  =  (a.  +  61)/(1  +  60),  i  =  1,2 


b^ '  M  0^/(1  +0^),  1  *  0,1,2  and  »  1 


For  the  deterministic  disturbance  we  can  find  such  a  polynomial 


P(z  )  =  (1  -  z  Yn  (l  -  2(  cos  2iri«)z  1  +  z  2),  j  >0 

i 


so  that 


P(z_1)Zt  -  0 


-1 . 


Thus,  multiplying  (A. 3)  by  P(z  )  gives  the  following  model 


-1 


-1 


-1 


A(z  ) (y  -  uy  +  6  )-  B(z  )(u  -  py  +  6  )  +  C(z  K 
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where 


A(z  *)  =  A'(z  Sp(z  S  =  1  +  ajZ  1  +  .  . .+  apz  ° 


B(  z  *)  =  B'(z  J)P(z  S  =  b  +  b  z  1  +  ...+  b  z  n 


C(z  )  =  C'(z  )B'(z  )P(z  )  *  1  +  CjZ  +  . . .+  cnz 

and  n  is  an  appropriate  positive  integer  which  represents  the  model 
order.  The  polynomials  of  A,  B  and  C  may  have  common  zeros  on  the 
unit  circle  because  of  the  periodic  disturbances. 

The  model  is  linear  only  if  overlap  factor  p  is  known 
analytically  or  experimentally.  A  special  case  is  the  plunge  cutting 
process  which  has  p  *  1. 


Appendix  2 


STC  With  Weighted  Minimum  Variance  Criterion 


The  control  criterion  of  weighted  minimum  variance 


2  2 

Min  E [ y  +  B'u  It],  B'  >  0 
t+k  t 


gives  the  following  control  law  for  (4.1) 


y  +  Bu  =  0,  B 

t+k|  t  t 


B'/g, 


The  self-tuning  algorithm  is  thus  (39J, 

(i)  choose  appropriate  B 

(ii)  estimate  the  model  parameters 
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*  forgetting  factor  (<1) 
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Appendix  3 


The  FORTRAN  Program  for  Digital  Simulation 


In  this  Appendix  we  give  the  FORTRAN  program  which  is  used  for 
the  digital  simulation  in  Chapters  5  and  6.  The  program  is  written 
for  the  proposed  GAC  system  with  multiprobe  measurement,  (4.33-35). 
By  modifying  the  program  so  that  the  probe  function 

pd<r‘>  - 1 

and  the  detuning  factor 


P  =*  1 


the  program  can  also  be  run  for  the  GAC  system  without  multiprobe 
measurement,  (4.11)  and  (4.12). 

Subroutine  SIML  is  the  simulated  plunge  cutting  process  as 

I 


described  in  Chapter  5.  The  adaption  mechanism  is  programmed 
according  to  the  minimum  variance  criterion  as  shown  in  Appendix  2. 
Subroutine  SQRTES  is  a  recursive  least  squares  routine  with  square 
root  factorization  [41] .  Three  types  of  the  varying  forgetting 

factors  are  programmed  in  subroutine  VFORG.  The  third  type,  which  is 
proposed  by  Fortescue  et  al.  [36],  needs 
FI  *  (constant  chosen  for  appropriate  adaption  speed) 

F2  *  1  + 


F3  =  e 


non 


PROGRAM  PLGMLP( INPUT, OUTPUT, PLOT, T AP E5 “INPUT, TAP E6 -OUTPUT) 

******** ************************************************* ********* 

ADAPTIVE  ROUNDNESS  CONTROL 

C  THIS  IS  A  TECHNIQUE  OF  ADAPTIVE  CONTROL  WITH  MULTIPROBE 
C  MEASUREMENT.  IT  HAS  TWO  MAJOR  PARTS: 

C  1.  RECURSIVE  PARAM  EST 

C  2.  IMPLICIT  MIN  VAR  CONTROL 

C  REF:  ASTROM  &  WITTENMARK  (1973) 

C  CLARKE  (1981) 

C  WHITEHOUSE(  1976) 

C  PH.D.  THESIS  BY  YHU-TIN  LIN,  AUG.,  1985,  PURDUE  UNIV. 

C  GEN  LIB:  IMSL 

C* ****************************************** *********************** 

c 

COMMON  RNS , FREQ, PHASE, AMPL, ALPH, BETA, GAMA ,JF ,JG 

COMMON  Y(IOOO) ,U( 1000) ,E( 1 000) ,X(20) , YS( 1 000) , FEED, F( 1 000) 

DIMENSION  A(10) ,B(10) ,C(10) ,YP(1000) 

DIMENSION  Z(20) ,THETA(20) ,THETB(20) , S( 1 20) , ZY( 20) , ZU( 20) ,ZD(20) 
C 

READ  200, IP, IPRINT,IFEED,K,N,NA,NEEDPLT,NT 
C 

Kl-K+l 

L-N+K+l 

LK-L+K 

LK1-L+K+1 

LN1-L+N+1 

Nl-N+1 

NONA+1 

NK-N+K 

NPAR»N+L+N 

NTO-NT-1 

NTK-NT+K 

NTP-NT~IP*2 

C 

C  *****  SETUP  OF  MULTIPROBE  CONFIG. 

C 

READ  200.JF.JG 
DP 1*2. *3. 1415926 

ALPH-SIN(DPI*(JG-JF)/IP)/SIN(DPI*JG/IP) 

BETA— 1.0 

GAMA“SIN(DPI*JF/IP)/SIN(DPI*JG/IP) 

C 

c  *****  SETUP  OF  SIMULATION  CONDITIONS 
C 

READ  210,  AMPL, FEED, PHASE, PSQR, RNS, STIFF 
FREQ*2.*(2.*3. 141 592654/ IP) 

READ  220, (A(J) ,J-1 ,NA) ,(B(J) , J-l ,NA) , ( C( J ) , J- 1 , NC) 

PRINT  230,  RNS, FREQ, AMPL, JF,JG,  ALPH, BETA, GAMA 
PRINT  240, NT, N, K, I P.NPAR.IFEED, FEED, STIFF, PSQR 
C 


o  non 


c  *****  SELECT  FORGETTING  FACTOR  AND  DETUNING  FACTOR 

JF0RG=>1  F1-F2/I 

2  Fl*F0RGETfF2 

3  FORTESCUE(  1981) 

READ  250, JFORG, FI , F2 , F3 , DTUNE 
PRINT  2 70, JFORG, FI, F2.F3, DTUNE 

C 

C  *****  INITIALIZATION  OF  PLOT 
C 

C  NEEDPLT>0  NO  PLOT 

C  1  TWO  POLAR  PLOTS 

C  2  ALL  PLOTS 

C  3  CONTROLLED  POLAR  &  TIME  SERIES  (UNSTABLE  SYSTEM) 

C  4  ONLY  TIME  SERIES 

C 

IF  (NEEDPLT.EQ.O)  GO  TO  30 
CALL  PLOTS 
CALL  FACTOR( .4) 

SFPOLA^40. 

XC-0. 

Y012.5 

YXL-4. 

CALL  PLOT( 12.5,0. ,3) 

CALL  PLOT( 12.5,16. ,3) 

IF  (NEEDPLT.EQ. I)  GO  TO  30 
C 

SFPARA-1. 

SFTIMD- 1  5. 

XINI-0. 

XLEW-10. 

YINI-4.2 

YLEN-2. 

THERO-2. 

THMIN--THERO/SFPARA 

CALL  LABEL(NT,XINI,XLEN,YINI,YLEN,THMIN) 

DT-XLEN/NT 

T-XINI 

TO-XINI 

UZERO-YINI+1. 

YZERO-YINI+3. 

YZERl>*YINI+5. 

DO  20  J-l.NPAR 
20  THETB(J)-0. 

C 

C  *****  TEST  UNCONTROLLED  SYSTEM 

C  IFEED  -  0,  NO  RETRACTION;  -1,  WITH  RETRACTION 

C 

30  CONTINUE 
DO  40  J-l.NA 
B(J)-B(J)*STIFF 
40  X(J)-0. 


o  n 


DO  50  1=1 ,  NTK 
YS( I)=FEED*I 

IF  (I.GT.NTP)  YS(I)=YS(NTP) 

IAP-I-IP 

IF  (IAP.LE.O)  GO  TO  50 
IF  (IFEED.EQ.O)  YSO-YS(IAP) 

IF  (IFEED.EQ.O)  YSO=FEED*IAP 
50  F ( I ) =Y S ( I ) -Y SO 
ICTL-0 

CALL  SIML( A,B,C,NT,NA,IP,  ICTL) 

IL-NT 

DO  60  1=1, NT 
F(IL)=PROBE(F,IL) 

IL-IL-1 

60  YP(I)=Y(I)-YS(I) 

NTR-NT-IP 

IF  (NEEDPLT.NE.l . AND.NEEDPLT.NE.2)  GO  TO  65 

NTR=NT-IP 

XC=XC+2. 5 

IFCOtt-O 

CALL  RONDPLT(YP,XC ,YC ,YXL,NTR,IP, RNS , IFCON, SFPOLA) 
65  IF  (NEEDPLT.LE.l )  GO  TO  70 

CALL  SPLOT(YP,NT,YZERU,TO ,DT, SFTIME) 

CALL  PL0T(XINI,THER0,3) 

C 

C  *****  INITIALIZATION  FOR  PARAM  EST 
70  PRINT  280 
ICTL*! 


READ  220, (THETA(J) ,J=1 ,NPAR) 

THETA(N1)=1 . 

DO  80  J=1 , NA 
80  X(J)=0. 

DO  90  1=2, LK1 

90  CALL  SIML(A,B,C,I,NA,IP, ICTL) 

CALL  INIT( YP, ZY, ZU, ZD , S ,N,K, I P, PSQR) 

*****  RECURSIVE  EST  BY  SQUARE  ROOT  METHOD  W.  FORGET  FACTOR 
C 

ZDK=ZU( 1 )-PROBE(U ,LK) 

DO  110  I=LK1, NTO 

CALL  VFORG(  FORGET, JFORG,  I ,  FI ,  F2  ,  F3  ) 

IY-I-K 

YP(I)=Y(IY)-YS(IY) 

CALL  MOVE ( ZY , K 1 , Z ,  1 ,N) 

CALL  M0VE(ZU,K1, Z,N1,L) 

CALL  M0VE(ZD,K1, Z,LN1,N) 

ZZ-ZDK 

IF  (K.NE.O)  ZZ=-ZD(K) 

PHI=PROBE(  YP,  I  )+ZZ 

CALL  SQRTES(PHI, Z, S, THETA, NPAR, ERR, FORGET) 
c  *****  DATA  UPDATE  &  CONTROL 
CALL  UPDAT( ZY, PHI,NK) 
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CALL  UPDAT(Z(J, 0.  ,LK) 

CALL  UPDAT(ZD, -ZDK,  NK) 

SUM-SCAPRD(  ZY, THETA, 1 , 1  ,N)+SCAPRD(ZU, THETA, 1  ,N1,L)+ 
&  SCAPRD( ZD, THETA, 1 ,LN1,N) 

1 1-1+1 
IK1-I1+K 
IKP-IK1-IP 
ZDK-F(IKl) 

IF  (IKP.GT.O)  ZDK-ZDK-PROBE(YP,IKP) 

ZU( 1)-(ZDK-SUM)/THETA(N1) 

U(I)-(ZU( 1 )-ZDK-PROBE(U, I ) ) /ALPH*DTUNE 
ZU(  1)-PR0BE(U,  I)+ZDK 


140 

**’■,« 

m 

150 

v*. 

200 

m*, 

210 

" 

220 

M 

rcr. 

230 

IYO-IY-1 

IF  (I.GT.IPRINT)  PRINT  260,  YP(  I) ,U( IYO) ,ERR,( THETA( J ) ,J-1 ,NPAR) 
IF  (NEEDPLT.LE.l)  GO  TO  100 

CALL  CPLOT( THETA, THETB,NPAR,T,DT,THERO,SFPARA) 

C  *****  COMPUTE  SYSTEM  OUTPUT 

100  CALL  SIML(A,B,C,I1 ,NA,IP, ICTL) 

110  CONTINUE 
C 

C  *****  PLOT  CONTROL  INPUT  &  OUTPUT 
C 

DO  120  I=NT,NTK 

IY-I-K 

IYO-IY-l 

YP(  I)=Y(IY)-YS(IY) 

120  PRINT  260,YP(I),U(IY0) 

IF  (NEEDPLT.EQ.0.0R.NEEDPLT.EQ.4)  GO  TO  130 

NTKP-NTK-IP 

XC-XC+5. 

IFCON-1 

CALL  RONDPLT( YP,XC ,YC ,YXL,NTKP, IP, RNS , IFCON, SFPOLA) 

130  IF  (NEEDPLT.LE.l)  GO  TO  150 
C 

DO  140  1*1 , NT 
IK-I+K 

YP(  I)-YP(  IK) 

CALL  SPLOT(YP,NT,YZERO,TO , DT, SFTIME) 

CALL  SPLOTCU, NT, UZERO, TO, DT, SFTIME) 

CONTINUE 

CALL  PL0T(0. ,0. ,999) 

STOP 


&  12,'  JG-' ,12, '  ALPtt-',F10.6,'  BETA-' ,F10.6,'  GAMA-' ,F10.6) 

240  FORMATC 1 X, 'NT- ' , 14 , '  N«' ,12, '  K-',I2,'  IP-', 12,'  NPAR  ',12, 

&  '  IFEEEH',12,'  FEE>',F5.4,'  STIFF-'  ,F5.2, '  PSQR-',F5.1) 

250  FORMATC 15, 4F1 0.0) 

260  FORMATC IX, 16F8. 5) 

270  FORMATC IX, 'TYPE  OF  FORGET:' , 12 , '  F1-',F8.5,'  F2-',F8.5,'  F3-', 


noon  o  o  o  o  on 


&  F8 . 5 ,  '  DETUNE*  ' ,  F8 . 5  ) 

280  FORMATC IX,  80( '-')/7 X, 'Y-' ,6X, 'U-' ,4X, 'ERR-' ,2X, 'THETA-' ) 

END 

SUBROUTINE  INIT(YP, ZY, ZU, ZD,S,N,K, IP.PSQR) 

Q  **************************************************************** 

C  INITIALIZATION  OF  VARIABLES 

*** ********* ******************** ******** A A A* A  4c  ** ** ************** 

COMMON  RNS, FREQ, PHASE, AMPL.ALPH, BETA, GAMA.JF.JG 

COMMON  Y( 1 000) ,U( 1 000) , E ( 1 000) ,X( 20) , YS( 1 000) , FEED, F( 1 000) 

DIMENSION  YP(IOOO) ,ZK(20) ,ZY(20) ,ZU(20) ,ZD(20) , S( 1 20) 

L-N+K+l 

NK-N+K 

Kl-K+1 

LK-L+K 

DO  2210  I-l.LK 
YP( I)=0. 

IK-I-K 

2210  IF  (IK.GT.O)  YP( I)-Y( IK)-FEED*IK 

LKK1-LK+K+1 
DO  2220  J-l.LK 
LKP-LKK1-IP 
ZK( J)-F(LKKl) 

IF  (LKP.GT.O)  ZK( J) -ZK( J)-PROBE( YP, LKP) 

2220  LKK1-LKK1-1 

***  INIT  Z( .)  &  X(  .) 

JY-LK 

DO  2230  J-l.NK 
JU-K1+J 

ZY(J)»PROBE(YP,JY)+ZK(JU) 

2230  JY-JY-1 

DO  2240  J-l.LK 
ZU(J)-PROBE(U,JU)+ZK(J) 

2240  JU-JU-1 

DO  2245  J-l.NK 
Jl-J+1 

2245  ZD(J)— ZK(J1) 

***  INIT  S( .)TRANSPOSE-  UPPER  TRIANGLE,  IN  ROW-WISE 
C 

NPAR-N+L+N 

JI-1 

DO  2250  J-l.NPAR 
DO  2250  I-J.NPAR 
S(JI)r 

IF  (I.EQ.J)  S( J I) -PSQR 
2250  JI-JI+1 
RETURN 


END 

SUBROUTINE  SIML(A,B,C,IT,N, IP, ICTL) 
c  **************************************************************** 

C  SIMULATED  PLUNGE  CUTTING  PROCESS 
C  ICTL-0  (UNCONTROL) 

C  ICTL-1  (CONTROL) 

*********** ******************* **************** ************ ****** 

C 

COMMON  RNS, FREQ, PHASE, AMPL.ALPH, BETA, GAMA, JF,JG 
COMMON  Y(IOOO) ,U( 1 000) , E( 1 OOO) ,X(20) ,YS(1000) ,FEED 
DIMENSION  A(20) ,B(20) ,C(20) 

11- IT 

12- IT 
Nl-N+1 

IF  (ICTL.NE.O)  GO  TO  2120 
11-2 

SEED-1 2347.D0 

CALL  GGNML( SEED, IT, E) 

RU-  .001 

SEED-127. DO 

CALL  GGNML(SEED.IT.U) 

DO  2110  I-l.IT 

U(I)-U(I)*RU 

IF  (I.GT.N1)  U(I)-0. 

2110  E(I)-E(I)*RNS 

Y(  1 ) -E  ( 1 ) +Y  S  ( 1 ) +AMP  L*  S I  N(  F  REQ4-P  H  ASE  ) 

C 

2120  CONTINUE 

DO  2150  1-11,12 
El-E(I) 

IN1-I-N1 

IF  (IN1.GT.0)  E1-E1+C(N1 ) *E( INI ) 

X(N1)=0. 

DO  2140  J-l , N 
IJ-I-J 

IF  (IJ.LE.O)  GO  TO  2130 
E1-E1+C(J)*E(IJ) 

Yl-Y(IJ) 

IJP-IJ-IP 

IF  (IJP.GT.O)  Yl-Yl-Y(IJP) 

X(N1)-X(N1)+B(J)*Y1 
2130  JN-N1-J 

X(N1)-X(N1)-A(JN)*X(J) 

Jl-J+1 

2140  X(J)-X(J1) 

Y( I )-E 1-X( N1 )+YS( I )+AMPL*S IN( FREQ*I+PHASE) 

IO-I-l 

Y(r)-Y(I)+U(IO) 

IBP-I-IP 

IF  (IBP.LE.O)  GO  TO  2150 
YB-Y( I )-Y( IBP) 

IF  (YB.LT.O.)  YB-O. 

Y( I)-YB+Y( IBP) 


2150  CONTINUE 
RETURN 
END 

SUBROUTINE  VFORG( FORGET,  JFORG, I , FI , F2 , F3) 

C  ******************************************************************* 

C  COMPUTE  TIME  VARYING  FORGET 

q  ******************************************************************* 

C 

IF  (JFORG. EQ.l)  FORGET* F1-F2/ I 
IF  (JFORG. EQ. 2)  FORGET=Fl*FORGETf F2 
IF  (JFORG. NE .3)  RETURN 
FORGET* l.-F3*F3/(Fl*F2) 

RETURN 

END 

SUBROUTINE  SQRTES(PHI, Z ,S , THETA, N,E , FORGET) 

C  *************************************************************** 

C  RECURSIVE  PARAM  ESTIMATION  OF  SISO  SYSTEM 

C  BY  SQUARE  ROOT  METHOD 

C  REF:  CLARKE  (1981) 

C 

C  *************************************************************** 

c 

DIMENSION  S(1 20) ,Z(20) ,GK(20) ,G(20) ,THETA( 20) 

C 

C  ***  CAL  PRED  ERROR 
C 

E=PHI-SCAPRD( Z, THETA,  1 ,1,N) 

C 

C  ***  CAL  SCALAR  DIVISOR-S IGSQ  AND  KALMAN  GAIN 
C 

SIGSQ-FORGET* FORGET 

K-N 

JI-1 

DO  3090  J-l.N 
3090  GK(J)=0.0 
C 

DO  3100  J=1,N 
F J-S CAPRD(  S  ,  Z  ,J  I ,  J ,  K ) 

K-K-l 

SIGSQ-SIGSQ^FJ*FJ 
DO  3100  W,N 
GK(I)-GK(  I)+FJ*S(JI) 

3100  JI-JI+1 
C 
C 

SIG-SQRT(SIGSQ) 

DO  3110  1=1, N 
3110  G(I)<K5K(  I)/SIG 
C 

C  ***  RECURSIVE  ALGORITHM  ON  S(.) 

C 

K-l 

DO  3130  J=1 , N 


A-SQRT(S(K)*S(K)-G(J)*G(J)) 

U»S(K)/A 

V— G(J)/A 

S(K)=A/FORGET 

IF(J.EQ.N)GO  TO  3130 

Kl-K+1 

Jl-J+1 

DO  3120  I=J 1 , N 
B-S(Kl) 

S(K1) “(U*S(K1)+V*G( I) ) /FORGET 
G(I)-U*G(I)+V*B 
3120  Kl-Kl+1 
K-Kl 

3130  CONTINUE 
C 

C  ***  CAL  NEW  PARAMETERS 
C 

A-E/SIGSQ 
DO  3140  1*1, N 

3140  THETA(I)-THETA(I)+GK(I)*A 
RETURN 
END 

FUNCTION  SCAPRD(  A1,A2,I1,I2,L) 

Q  ★***★**★★*★ ********<k  ********************************** ********* 
C  CAL.  SCALAR  PRODUCT  OF  TWO  COLUMN  MATRICES 

C  *************************************************************** 

C 

DIMENSION  Al(l 00) ,A2(100) 

SCAPRI>0.0 

Jl-Il 

J2=*I2 

DO  2410  1*1 ,L 

S CAPRD-S CAPRDf  A1 ( J 1 ) *A2 ( J 2) 

J  1-J1+1 
2410  J2*J2+1 
RETURN 
END 

SUBROUTINE  MOVE(A,IA,B,IB,LEN) 

Q  ********************************************************** 

C  MOVE  A  TO  B 

q  ********************************************************* 

C 

DIMENSION  A(50) ,B(50) 

JA-IA 

JB-IB 

DO  2510  I-l.LEN 
B(JB)-A(JA) 

JA-JA+1 
2510  JB-JB+1 
RETURN 
END 

SUBROUTINE  UPDAT( Z ,Y, NPAR) 

C  *************************************************************** 


o  o  n  o  non 


C  SHIFT  I/O  DATA  (FIFO) 

C  *************************************************************** 

c 

DIMENSION  Z(20) 

J-NPAR 

DO  2310  I-2.NPAR 
J 1-J-l 
Z(J)-Z(J1) 

2310  J-Jl 
Z(  1 )  =*Y 
RETURN 
END 

FUNCTION  PROBE(Y.I) 

C  *************************************************************** 

FILTERED  OUTPUT  OF  THREE  PROBES 

*************************************************************** 


COMMON  RNS, FREQ, PHASE, AMPL, A, B,C,JF ,JG 
DIMENSION  Y ( 1 000) 

IJF-I-JF 
I JOI-JG 
PROBE-Y(I)*A 

IF  (IJF.GT.O)  PROBE*>PROBE+Y( IJF)*B 
IF  (IJG.GT.O)  PROBE>PROBE+Y(IJG)*C 
RETURN 
END 

S UB ROUTINE  ACCUR(Y, I , N , D IMERR .OUTRND, FINISH) 

****************************************************************** 

EVAL.  ACCURACY  OF  MACHINED  PART 

****************************************************************** 


DIMENSION  Y(1000) 

SUM-0. 

SUMSQ-0. 

YMAX-0. 

YMIN-O. 

Kl-I+1 

K2-I+N 

DO  2520  K-K1.K2 

IF  (Y(K) .GT.YMAX)  YMAX-Y(K) 

IF  (Y(K) .LT.YMIN)  YMIN-Y(K) 

SUM-SUM+YOO 

2520  SUMSQ-SUMSQ*-Y(K)*Y(K) 

DIMERR-SUM/N 

FINIS&-SQRT(SUMSQ/N-DIMERR*DIMERR) 

OUTRNI>YMAX-YMIN 

PRINT  *,'OUTRND' .OUTRND,'  FINISH' .FINISH,'  DIMERR' .DIMERR 

RETURN 

END 

SUBROUTINE  CP LOT( THETA, THETB ,NPAR,T,DT, OFFSET, SF, 

C  *************************************************************** 

C  PLOT  OF  PARAMETERS 

r  *************************************************************** 


o  o  o 


DIMENSION  THETA(20) ,THETB(20) 

T1*T 

T*T+DT 

THETB(  1 )  =THETA(  1 ) 

DO  4110  J=1 , NPAR 

CALL  P LOT( T , OFFSET+T HETA( J  )  *S  F  ,2  ) 

J1*J+1 

IF(Jl.LE.NPAR)  GO  TO  4100 

Jl*l 

Tl-T 

4100  CALL  PLOT( T 1 ,OFFSET+THETB( J 1) *SF ,3) 

4110  THETB( J  1)=THETA(J  1) 

RETURN 

END 

SUBROUTINE  SPLOT(Y,NT, ZERO, TO , DT, SF) 

Q  *****************************  **************************  ******* 

PLOT  OF  SINGLE  TIME  SERIES 

**************************** ********************************** 


DIMENSION  Y(2000) 

CALL  PLOT( TO, ZERO, 3) 

T-TO 

DO  4115  1*1, NT 

T-T+DT 

Y(I)*Y(I)*SF 

IF(ABS(Y(I)).GT.1,0)Y(I)”S IGN( 1 . ,Y( I) ) 

4115  CALL  PLOT(T,Y(I)+ZERO,2) 

RETURN 

END 

SUBROUTINE  RONDPLTCY.XC ,YC ,YXL, I ,N, RNS ,IFCON, SF) 

C  ************************************************************* 
C  POLAR  PLOT  OF  FORM  GEOM.  AND  EVAL.  OUT-OF-ROUNDNESS 
C  ************************************************************* 

C 

DIMENSION  Y(1000) 

CALL  AX I SX( XC , Y C-Y XL/ 2 , '  ' ,1 ,YXL,90. ,0. ,  .1 ,  .1,31) 

CALL  AX  I SX( XC- Y  XL/ 2 , YC  , '  ',1 ,YXL, 0.0,0. , .1 , .1 ,31) 

R*l .  5 

DANG” 2. *3. 141592 6/ N 
CALL  PLOT(XC,YC,-3) 

CALL  CIRCLE( R,DANG,N) 

CALL  NEWPEN( 2 ) 

C 

A”Y ( I ) *SF+R 
CALL  PL0T(A,0.  ,3) 

11*1+1 

I2-I+N-1 

DO  3500  K-11,12 

R1*R+Y(K)*SF 

ANG”DANG*(K-I) 

X  1*R1*C0S(ANG) 

Y 1*R1 *S IN( ANG) 


3500  CALL  PLOT(X 1 , Y 1 , 2 ) 

CALL  PLOT(A,0. ,2) 

CALL  NEWPEN(l) 

CALL  ACCUR(Y, I, N,DIMERR,OUTRND, FINISH) 

HIGH*. 2 
XS=-2. 

XN-.5 

YSN-3.0 

IF  (IFCON.EQ.O)  CALL  SYMB0L(XS ,YSN, HIGH,'  (UNCONTROLLED) ' ,0. , 1 7) 
IF  (IFCON.EQ.l)  CALL  SYMB0L(XS .YSN.HIGH, '  (CONTROLLED) ' ,0. , 1 7) 
YSN-YSN-HIGH*1.5 

CALL  SYMBOL(XS, YSN.HIGH, 'DIMENS  ERROR:' ,0. ,1 2) 

CALL  NUMBER(XN,YSN,HIGH,DIMERR,0. ,5) 

YSN“YSN-HIGH*1 . 5 

CALL  SYMBOL(XS, YSN.HIGH, 'OUT-OF-ROUND:' ,0. ,12) 

CALL  NUMBER(XN,YSN,HIGH,OUTRND,0. ,5) 

IF  (RNS.EQ.O)  GO  TO  3550 
YSN-YSN-HIGH*1.5 

CALL  SYMBOL(XS, YSN.HIGH, 'FINISH  (RMS) : ' ,0 . , 1 2) 

CALL  NUMBER(XN, YSN.HIGH, FINISH, 0. ,5) 

3550  CALL  PLOT(0. ,0. ,3) 

CALL  PLOT(-XC ,-YC ,-3) 

RETURN 

END 

SUBROUTINE  C IRCLE( R.DANG, N ) 

CALL  PLOT( R ,0 . ,3) 

DO  3600  1=1, N 

ANG-DANG*I 

X=R*COS(ANG) 

Y=R*S IN( ANG) 

3600  CALL  PLOT(X,Y,2) 

RETURN 

END 

SUBROUTINE  LABEL( NT.XIO , XLO , Y 10 , YLO , THMIN) 

c  *************************  **  *******************  ********  it  Irk  it  ft  irk* 

C  PLOT  OF  LABELING  AND  COORDINATE  FRAME 

C  **************************** ********** * A********************** 

C 

DX-NT/XLO 

XI-XIO 

XL-XLO 

YI-YIO 

YL-YLO 

Y=*Y  I 

CALL  AXIS(XI,YI, '  ' ,-l , XL, 0. ,0. ,DX, 20) 

DO  4220  J-1,3 

CALL  AXIS(XL,Y,'  ' , 1 , YL, 90. ,0 . , 1 . , 1 6) 

CALL  AXIS(XI,Y, '  ' ,-l , YL, 90. ,0. , 1 . , 1 6) 

DO  4220  K-1,2 
Y-Y+l. 

CALL  PLOT(XI , Y , 3) 

CALL  PLOT(XL,Y,2) 

4220  CONTINUE 
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YI-O. 

YL-4. 

DY=-THMIN*2./YL 
CALL  PL0T(XI,YL,3) 

CALL  PL0T(XL,YL,2) 

CALL  AXIS(XI , Y I, '  ' ,-l , XL, 0. ,0. ,DX, 4 ) 

CALL  AXIS(XL, YI, '  ' , 1 , YL, 90. ,0. , 1 . , 1 6) 

CALL  AXIS(XI,YI, 'THETA', 5,YL, 90. ,THMIN,DY, 7 ) 

XS-.50 

YS-Y-.4 


X  N“  -  .  3 
YN-Y-1. 

HIGH-. 2 

CALL  SYMBOL(XS,YS, HIGH, 'UNCONTROL  RESPONSE' ,0 ., 1 8) 
CALL  NUMBER(XN,YN,HIGH,0. ,0. ,0) 

YS-YS-2. 

YN-YN-2. 

CALL  SYMBOL(XS,YS, HIGH, 'CONTROL  RESPONSE' ,0 ., 1 6) 
CALL  NUMBER(XN,YN,HIGH,0. ,0. ,0) 

YS-YS-3. 

YN-YN-2. 

CALL  SYMBOL(XS,YS, HIGH, 'CONTROL  INPUr,0.,13) 

CALL  NUMBER(XN,YN,HIGH,0.  ,0.  ,0) 

RETURN 

END 
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C.  R.  Liu,  Professor  of  Industrial 

Engineering .  Principal  Investigator 

K.  S.  Fu,  Goss  Distinguished  Professor 

of  Engineering  (Elec. Eng.) .  Faculty  Associate 

J.  Modrey,  Professor  of  Mechanical 

Engineering .  Co— Principal  Investigator 


A.  L.  Sweet,  Professor  of  Industrial 
Engineering . 


W.  Stevenson,  Professor  of  Mechanical 
Engineering . 


J.  J.  Talavage,  Professor  of  Industrial 
Engineering . 


Co-Principal  Investigator 


Faculty  Associate 


Faculty  Associate 


R.  Hannam",  Visiting  Associate  Professor 

of  Industrial  Engineering .  Faculty  Associate 

W.  Johnson* ",  Visiting  Professor  of  Industrial 

Engineering .  Faculty  Associate 

A.  Shumsherrudint,  Visiting  Associate  Professor 

of  Industrial  Engineering .  Faculty  Associate 
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